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ABSTRACT 

We present the mini-proceedings of the workshop on "Parton fragmentation in the vacuum and in the 
medium" held at the European Centre for Theoretical Studies in Nuclear Physics and Related Areas 
(ECT*, Trento) in February 2008. The workshop gathered both theorists and experimentalists to discuss 
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the current status of investigations of quark and gluon fragmentation into hadrons at different accelerator 
facilities (LEP, B-factories, JLab, HERA, RHIC, and Tevatron) as well as preparations for extension of 
these studies at the LHC. The main physics topics covered were: (i) light-quark and gluon fragmentation 
in the vacuum including theoretical (global fits analyses and MLLA) and experimental (data from e^e^ , 
p-p, e-p collisions) aspects, (ii) strange and heavy-quark fragmentation, (iii) parton fragmentation in 
cold QCD matter (nuclear DIS), and (iv) medium-modified fragmentation in hot and dense QCD matter 
(high-energy nucleus-nucleus collisions). These mini-proceedings consist of an introduction and short 
summaries of the talks presented at the meeting. 
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1 Introduction 



The transition from coloured quarks and gluons to colourless hadrons - the so-called fragmentation or 
hadronization process - is a Quantum Chromodynamics (QCD) phenomenon with many important the- 
oretical and phenomenological implications for the physics at high-energy colliders. At large hadron 
fractional momenta z = /^hadron/T^parton ^ O-Ij the fragmentation functions (EEs) of partons into hadrons 
or into photons obey DGLAP evolution equations and are obtained from global analyses (fits) of var- 
ious experimental hadron and photon production data. The EEs are, for example, a basic ingredient 
for the calculation of the production of high transverse-momentum particles at collider energies within 
perturbative QCD. At small z, successful QCD resummation techniques (e.g. the Modified Leading Log- 
arithmic Approximation, MLLA) have been developed to understand and describe the evolution of a 
highly-virtual time-like partons into final hadrons. Most of the available data used in the study of the 
fragmentation process comes from qq{g) production at e^e^ colliders (LEP, SLC), although data from 
deep-inelastic e-p scattering (DIS) at HERA has also been used especially in the heavy-quark sector 
New high-precision flavour-identified hadron data from the B-factories (BELLE, BABAR), from DIS 
(HERA), and from hadronic colliders (RHIC, Tevatron), can further help to constrain the EEs and, in 
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fine, to improve the perturbative calculation of hadron and photon production at the LHC (in a way simi- 
lar to the recent developments in improved global-analyses of parton distribution functions). In the near 
future, the vast kinematical range opened in p-p collisions at LHC energies and luminosities will open 
up new channels for the study of fragmentation functions. 

The study of possible modifications of the parton fragmentation processes in heavy-ion collisions is also 
an important tool for the determination of the thermodynamical and transport properties of the dense 
QCD matter produced at RHIC and LHC energies. Since the discovery of the suppression of high-pj- 
hadrons in central Au-Au collisions at RHIC, a lot of effort has been devoted to the understanding of the 
propagation of partons in QCD media (cold nuclear matter or quark-gluon-plasma) as well as to under- 
stand such a "medium-modified" fragmentation mechanism. These phenomenological studies - which 
often use techniques originally developed "in the vacuum" (e.g. MLLA) - are well supplemented by a 
wealth of new data from DIS on nuclear targets (HERMES at HERA, and CLAS at JLab) and in heavy- 
ion collisions (RHIC). 

A few months before the start of the LHC, it seemed a timely moment to have a workshop, gath- 
ering both theorists and experimentalists, to discuss the current status of investigations concerning the 
fragmentation mechanisms of gluons, light-quarks and heavy-quarks into hadrons and photons, as well 
as to discuss their imphcations for the upcoming p-p and A-A experimental programmes at the LHC. A 
meeting was organized at the ECT* (Trento) from February 25-29, 2008, with partial financial support 
from the center. The meeting had 50 participants whose names and institutes are listed in Appendix A. 
There were 43 presentations of various lengths. Ample time was left for discussions after each talk. The 
talks and discussions were organized around the following main topics: 

• Light-quark and gluon fragmentation in the vacuum (e~^e~, p-p, e-p collisions): global fits analy- 
ses and MLLA. 

• Strange and heavy-quark fragmentation. 

• Parton fragmentation in cold QCD matter. 

• Medium-modified fragmentation in hot and dense QCD matter. 

This mini-proceedings include a short summary of each talk including relevant references, the list of 
participants and the workshop programme. We felt that such a format was more appropriate a format 
than full-fledged proceedings. Most results are or will soon be published and available on arXiv. Most 
of the talks can also be downloaded from the workshop website: 

http : / / arleo . web . cern . ch/ arleo/f f _vacuum_medium_ect08 

We thank the ECT* management and secretariat, in particular Cristina Costa, for the excellent organiza- 
tion of the workshop and all participants for their valuable contributions. We believe that this was only 
the first workshop of this kind and look forward to similar meetings in the future. 



FRANfois Arleo, David d'Enterria 



3 



2 Theoretical aspects of parton fragmentation in the vacuum: global fits 
and MLLA 



Vacuum fragmentation functions: a common interface? 

Frangois Arleo, David d'Enterria 

A lot of effort has been invested over the past few years in order to gather all existing sets of Par- 
ton Distributions Functions (PDFs) under a common interface. This idea was developed in the Les 
Houches Workshop 2001 and later gave birth to the LHAPDF interface lUlIll, aimed to be the successor 
of PDFLIB |31. The essential points of LHAPDF, summarized in [111, are: 

• PDF sets including their uncertainties should be easily handled, such that any user can quantify 
the influence of PDF uncertainties on a given observable; 

• QCD evolution codes can be compared easily; 

• PDF sets are given under a parametrized form at a soft scale Qq and not through heavy grids (yet 
this is now possible in practice for older sets); 

• LHAPDF is modular and therefore flexible enough to facilitate the inclusion of new evolution 
codes or new FF sets. 

As shown in the workshop, there is also an important ongoing activity in the FF parametrizations 
community with e.g. the AKK08, DSS, and HKNS sets recently released, which complement and refine 
previously existing analyses (AKK, BFG, BKK, KKP, Kretzer, . . . ). All the existing sets with their main 
characteristics are summarized in Tabled 



Name 


Ref. 


Species 


Error 


Zmin 




(GeV2) 


AKK 


H 




no 


0.1 


2- 


-4-104 


AKK08 





■K^,K^,K^, p,pA,A 


yes 


0.05 


2- 


-4-104 


BKK 


El 


n++Ti-, 71°, K++K-,K^ + K^, h+ + 
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0.05 


2 


-200 


BFG 
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no 


10-3 


2- 


1.2 -104 


BFGW 


m 




yes_^ 


10-3 


2- 


1.2 -104 


CGRW 


m 




no 


10-3 


2- 


1.2 -104 


DSS 


nom 


71''', K^, p, p, h'^ 


yes_^ 


0.05-0.1 


1 


- 10^ 


DSV 


113 


polarized and unpolarized A 


no 


0.05 


1 


-104 


GRV 


El 


Y 


no 


0.05 




> 1 


HKNS 


HI 


Ti^, Ti^, K^, K'^^ + k^,n,p + p 


yes 


0.01 - 1 


1 


-108 


KKP 




71+ + 7I-, 7i", K+ +K-, +K^,p + p,n + n, h+ + h- 


no 


0.1 


1 


-104 


Kretzer 


m 


Ti^,K^,h+ + h- 


no 


0.01 


O.f 


?-10^ 



Table 1 : Main characteristics of fragmentation function (FF) sets obtained from global fit analyses. 



Errors on the parameters and their correlation are given. 

^Typical uncertainties for "truncated" energy fractions Jo 2z£'f (z, Q^)dz studied for all flavors and hadron species using the 
Lagrange multiplier method. 
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Just like parton densities, the first error analyses in the fragmentation functions have also just ap- 
peared. In this context, it appears sensible to examine whether an effort similar to LHAPDF should be 
performed for the different FF sets. Ideally, a common interface would facilitate the systematic com- 
parison between the different sets, in order to explore for instance the theoretical uncertainties of NLO 
calculations for hadron production in p-p collisions at the LHC associated with the FFs (which are usu- 
ally much larger than the uncertainties due to the PDFs). As a very first step, a simple wrapper program 
should be developed to call any parton-to-hadron FF for any given species. Also, the latest produced FF 
parametrizations should be made available in a common public web-page such as the existing one from 
Marco Radici ifTTl . 

Albino-Kniehl-Kramer FFs: Improvements from new theoretical input and experimental data 

Simon Albino 

We present results from an update (AKK08) f5l of the previous AKK fragmentation function (FF) 
sets for charged pions, charged kaons, (anti)protons, neutral kaons and (anti)lambdas with additional 
theoretical and experimental input. We incorporate hadron mass effects 118,1 . and fit the hadron mass in 
the case of the e^e^ calculation which also has the effects of subtracting out other low ^/s and small-.x 
effects beyond the fixed order approach, such as higher twist, small-^; logarithms, etc. In the case of 
the baryons, the fitted masses are about 1% above the true masses, which is consistent with scenarios 
in which the baryons are produced mainly from direct partonic fragmentation with a small contribution 
from decays from slightly heavier resonances. A greater excess is found for the pion mass, suggest- 
ing contributions to the sample from decays of heavier particles such as p(770). Charged and neutral 
kaon masses are significantly below their true masses. A possible explanation for this is that there are 
significant contributions from complicated decay channels, such that the direct partonic fragmentation 
approach is insufficient. For this reason we do not impose SU(2) isospin symmetry of u and d quarks 
between charged and neutral kaons, which we do for pions. We implement large-^; resummation in 
the quark coefficient function of e^e^ reactions using the results from Ref. L19], since this is a simple 
improvement which modifies the cross section over the whole range in x that we constrain. Large-x 
resummation is also implemented in the DGLAP evolution of the FFs ||5l. This results in a significant 
improvement in the fit for charged kaons, (anti)protons and (anti)lambdas, while is essentially un- 
changed for neutral kaons and charged pions. 

In addition to the constraints of the previous AKK fit, we impose further constraints on the charge- 
sign unidentified FFs from the data for single inclusive production of identified particles {pj > 2GeV/c) 
from RHIC Il20ll2 1 112211231124 ] . the Tevatron [25] and electron-positron reactions below the Z pole mass 
and in the range 0.05 < x < 0.1. While the untagged measurements from electron-positron reactions 
provide excellent constraints for the sums of the charge-sign unidentified FFs for quarks of the same 
electroweak charges, they do not constrain the remaining degrees of freedom at all. As in the previous 
AKK fit, these were constrained using quark tagged data, while in the new AKK fit additional constraints 
are provided by the data from RHIC. These data are also much more sensitive to gluon fragmentation 
and impose (exclusively) new constraints on the charge-sign asymmetry FFs. Normalization errors were 
treated as systematic effects, i.e. these errors were incorporated via a correlation matrix. Their weights 
were fitted analytically and independently of the fit in order to further ascertain the quality of the fit, and 
their magnitudes were typically found to lie in the reasonable range of - 2. 

While the results for the fitted masses suggest that the baryons are the best candidates for studying 
direct partonic fragmentation, there unfortunately exist some inconsistencies between the calculation and 
the measurements of the inclusive production of these particles at RHIC: The description of the STAR 
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data for A/A fails, while the contribution from the initial protons' valence d quarks to the charge-sign 
asymmetry for p j'p from STAR is negative. Furthermore, while the contributions to the production from 
the fragmentation of the valence and sea quarks of the initial protons for the data from RHIC exhib- 
ited the expected behaviour, the contribution from the valence d quark fragmentation to the charge-sign 
asymmetry in (anti)proton production is negative. All these issues would be better understood in the 
context of an error analysis of the FFs. 

We compared our sets to the recent DSS OOl / DSV Ull and HKNS HH ones, and typically found 
reasonable agreement for favoured FFs but not for unfavoured, and large discrepancies exist at large-^; in 
some cases. Future hadron identified data from BABAR, CLEO, HERA and RHIC will help to clarify 
these issues and significantly reduce the large uncertainties in much of the FF degrees of freedom. 

De Florian-Sassot-Stratmann (DSS) global QCD analysis of fragmentation functions 

Marco Stratmann 

We present new sets fragmentation functions for pions, kaons lITOl . protons, and unidentified charged 
hadrons ifTTTl obtained in NLO combined analyses of single-inclusive hadron production in electron- 
positron annihilation, proton-proton collisions, and deep-inelastic lepton-proton scattering. At variance 
with previous fits, the present analyses take into account data where hadrons of different electrical charge 
are identified, which allow to discriminate quark from anti-quark fragmentation functions without the 
need of non-trivial flavor symmetry assumptions. The resulting sets are in good agreement with all data 
analyzed. The success of the global analysis performed here, including observables other than e^e^ 
annihilation, stands for an explicit check of factorization, universality, and the underlying framework of 
perturbative QCD. 

Increasingly accessible hadron production data from proton-proton collisions and hadron multiplic- 
ities coming from semi-inclusive DIS should not be disregarded as they offer a crucial piece of comple- 
mentary information that reduces significantly the uncertainties of the resulting fragmentation functions. 
For instance, we find that the new DIS pion and kaon multiplicities provided by the HERMES experi- 
ment effectively constrain the separation between favored and unfavored distributions, a separation that 
was either not implemented in previous sets or it was based on certain assumptions. The most recent 
RHIC results provide stringent constraints on the gluon fragmentation function and, in general, on the 
large z behavior of the other distributions. 

The implementation of the minimization in our global analysis is numerically fast and efficient 
and can be straightforwardly expanded to any future set of hadron production data. With the help of the 
Mellin moment technique ll26l . the entire analysis was consistently performed at NLO accuracy without 
resorting to often used approximations for NLO hard scattering cross sections. For completeness we 
also provide LO sets. An extensive use of the Lagrange multiplier technique [27] is made in order to 
assess the typical uncertainties in the extraction of the fragmentation functions and the synergy from the 
complementary data sets in our global analysis. 

Hirai-Kumano-Nagai-Sudoh (HKNS) fragmentation functions and proposal for exotic-hadron search 

Shunzo Kumano 

Fragmentation functions and their uncertainties are determined for pions, kaons, and protons by a 
global y} analysis of charged-hadron production data in electron-positron annihilation and by the Hes- 
sian method for error estimation 1,14] . The results indicate that the fragmentation functions, especially 
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gluon and light-quark fragmentation functions, have large uncertainties at small . There are large 
differences between widely-used functions by KKP (Kniehl, Kramer, and Potter) [151, AKK (Albino, 
Kniehl, and Kramer) [4|, and Kretzer [16|; however, they are compatible with each other and also with 
our functions if the uncertainties are taken into account. We find that determination of the fragmentation 
functions is improved in next-to-leading-order (NLO) analyses for the pion and kaon in comparison with 
leading-order ones. Such a NLO improvement is not obvious in the proton. Since the uncertainties are 
large at small (Jp-, the uncertainty estimation is very important for analyzing hadron-production data at 
small Q^- or pj {Q^ , pj » M|) in lepton scattering and hadron-hadron collisions. A code is available for 
general users for calculating the obtained fragmentation functions llT4l . 

Next, it is proposed that fragmentation functions should be used to identify exotic hadrons (2E\. 
As an example, fragmentation functions of the scalar meson /o(980) are investigated. The /o meson is 
considered as a candidate for an exotic hadron beyond the usual 1717 configuration because its strong-decay 
width is much larger than the experimental one |[29l if it is an ordinary qq meson. The radiative decay 
(|) —>■ /oY was proposed to find its internal structure f3U\, and subsequent measurements indicated tetra- 
quark or KK structure. However, its structure has not been clearly determined yet. Here, we investigate 
the possibility of ordinary qq, ss, tetra-quark, KK molecule, or glueball via the fragmentation functions 
of /o. It is pointed out that the second moments and functional forms of the u- and 5-quark fragmentation 
functions can distinguish the tetra-quark structure from qq. By the global analysis of /o(980) production 
data in electron-positron annihilation, its fragmentation functions and their uncertainties are determined. 
It is found that the current available data are not sufficient to determine its internal structure, while precise 
data in future should be able to identify exotic quark configurations. 

Main sources of uncertainty in quark and gluon fragmentation functions into hadrons and photons 

Monique Werlen 

A precise knowledge of the QCD hadron and photon production is necessary while searching for 
e.g. a low mass Higgs signal at LHC and Tevatron. Uncertainties on this background must therefore be 
asserted with care, in particular those due to the measurement of fragmentation functions. In the calcu- 
lation of inclusive hadron and y cross sections in pp collisions with the PHOX package - a set of NLO 
parton-level event generators for large pj PHOton, hadron and/or jet X-sections 1 3 1 1 - uncertainties in 
fragmentation functions (FFs) do matter compared to those from scales and parton distribution functions 
(PDF) at RHIC and LHC energies ll32l . FFs are obtained from fits to inclusive cross sections assuming a 
functional form at a given scale MfQ and evolving it to the scale of the data. The sources of uncertainties 
are both experimental and theoretical. Experimental uncertainties are due to statistical errors, in partic- 
ular at large momentum fraction z of the hadron in e^e^ data. But also systematic effects contribute 
such as in the data normalization or in the extrapolation (evolution) needed to cover the full z (Q^) range. 
Matching theory and data due to binning and cuts settings may also be an issue. Theoretical uncertainties 
come from the choice of the functional form for the selected z range, the choice of scales, the order of 
the theory (leading, next- to-leading NLO, resummation) and further parameters like PDFs and a^. 

As shown for the first time in [Ol, NLO fits to e^e^ — > n^X data from PEP and PETRA constrain the 
quark fragmentation into neutral pions while those to pp ^ n^X from ISR and UA2 constrain the gluon 
fragmentation. BFGW [ 8 1 sets of FF into unidentified charged hadrons have been obtained from LEP 
and PETRA data choosing optimized scales: a rough approximation of the large-z resummation ( |[T9l ). 
Gluon parameters are quite sensitive to the functional form chosen for the quark distribution. A full 
statistical error analysis is performed while estimation of the "theoretical error" as previously detailed 
can be obtained by comparison to other FF sets. Important discrepancies with the BKK set L6J have been 
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shown, specially at large z. while the gluon FFs at z > 0.5 is found much higher, in accordance with the 
UAl data, than in the Kretzer set fTF]. Large scale instabilities (specially at low ^/s and low pr) affect 
the phenomenology of inclusive production of pion [33]. Theoretical estimates rely on extrapolations of 
the FFs outside (0.75 <x < 0.9) of the region where they are actually constrained (0.1 < z < 0.7) by the 
data ll34l . A prehminary study |[35l with JETPHOX shows that the pj imbalance hadron-jet correlation 
cross-sections at RHIC energy may be used to constrain the FFs into hadrons in the high z region. For 
prih) > 25 GeV/c and prijet) > 30 GeV/c, the three set (BFGW, KKP, Kretzer) give quite different 
predictions: BFGW is systematically higher than KKP while Kretzer is a factor two lower. 

Thanks to recent data from PHENIX and DO spanning a large xj = 2pT / y/s range, inclusive prompt 
photons cross sections from ^/s = 23 GeV to = 1 .96 TeV are now well understood in the NLO QCD 
framework [36 1. Photons can be produced either directly or via a parton-to-photon FF. The latter be- 
comes important at low pj and high ^/s. Uncertainties on the parton-to-photon FFs are mainly due to the 
uncertainty on the gluon FF as quantified by the BFG sets I and II [7]. BFG-I reduces the cross section 
by up to 10% (a factor 2.5) at /^r = 3 GeV at RHIC (LHC) energy compared to BFG II [32]. As shown 
in |[35l the photon-jet correlation at RHIC energy may also be used to constrain the photon FF. Fixing 
the jet and varying the photon momenta allows for a direct measurement of the photon fragmentation at 
low z, a region barely accessible in the LEP experiments. 

In the calculation of inclusive hadron and photon cross sections in pp collisions with the PHOX 
programs - NLO event generators (parton level) for large pt PHOton (hadron or jet) cross-sections (X- 
sections) ifSTl -, uncertainties on fragmentation functions (FF) do matter compared to those from scales 
and parton distribution functions (PDF) at RHIC and LHC energies [32]. There are various sources 
of uncertainties in the determination of FF from fits to inclusive cross sections where one assumes a 
functional form at a scale Mpo and evolves it to the scale of the data. Experimental uncertainties are 
due to statistical errors (FF with a large momentum fraction z of the parton in e^e^ data), systematic, 
normalization, the probed z and range as well as the sensitivity to quark and gluon FF. Matching 
theory and data (binning, cuts) may also be an issue. Theoretical uncertainties come from the choice of 
the functional form, the z range of assumptions, the choice of scales, the order of the theory (leading, 
next-to-leading NLO, resummation) and further parameters like PDFs and tty. 

NLO fits to e+e^ n^X data from PEP and PETRA constrain the quark fragmentation into neutral 
pions while NLO fits to pp n^X from ISR and UA2 constrain the gluon fragmentation BFGW [HI 
sets of FF into unidentified charged hadrons have been obtained from LEP and PETRA data choosing 
optimized scales (a rough approximation of the large z resummation performed later in [|19.1 ). Gluon 
parameters were quite sensitive to the functional form chosen for the quark distribution. A full statistical 
error analysis is performed while estimation of the "theoretical error" embedded in the parametrizations 
and which come from various theoretical assumptions can be obtained by comparing with different sets: 
There are important discrepancies, specially at large z, with the BKK set [6] while the very low gluon 
FF at z > 0.5 from the Kretzer set [16] seems unfavored by UAl data. Large scale instabilities (spe- 
cially at low ^/s and low pr) affect the phenomenology of inclusive production of pion ll33l . Theoretical 
estimates rely on extrapolations of the FF outside (0.75 — 0.9) of the region where they are actually con- 
strained by the data (0.1 < z < 0.7) from which they are extracted (34]. A preliminary study [35] with 
JETPHOX shows that hadron-jet correlation at RHIC energy may be used to constrain the FF into hadron 
in the high-z region. With prih) > 25 GeV/c and prijet) > 30 GeV/c, BFGW is systematically higher 
than KKP lITSi while Kretzer is almost a factor two lower than BFGW. 

Thanks to recent data from PHENIX and DO spanning a large xj = lpT j range, inclusive prompt 
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photons cross sections from ^/s = 23 GeV to = 1 .96 TeV are now well understood in the NLO QCD 
framework 1*361. Photons can be produced either directly or via a parton-to-photon FF. The later become 
important at low px and high ^/s. Uncertainties on the photon FF are mainly due to the uncertainty on 
the gluon FF as quantified by the BFG sets I and II |7]. BFG I reduce the cross section by up to 10% (a 
factor 2.5) at pr = 3 GeV/c at RHIC (LHC) energy compared to BFG II (321. A preUminary study ISU 
with JETPHOX shows that photon-jet correlation at RHIC energy may be used to constrain the photon 
FF. With the momentum of the jet fixed, varying the momentum of the non isolated photon allows for a 
direct access to the photon fragmentation at low z, a region barely accessible in the LEP experiments. 

Time-like splitting functions at NNLO in QCD 

Sven-Olaf Moch 

We have employed relations between space-like and time-like deep-inelastic processes in perturbative 
QCD to calculate the next-to-next-to-leading order (NNLO) contributions to the time-like non-singlet 
and singlet quark-quark and gluon-gluon splitting functions for the scale dependence (evolution) of the 
parton fragmentation distributions Dj {x,Q^). The evolution equations read 

= /^f Pj(z,a.(G^))Dt(^,e^), (1) 

where x denotes the fraction of the momentum of the final-state parton / carried by the outgoing hadron 
h, is a time-like hard scale and summation over j = q, g is understood. The time-like splitting 
functions governing Eq. ([T]) admit an expansion in powers of the strong coupling tts, 

Pji[x,a,{Q')) = a,pf\x)+alpf\x)+alpf\x) + ... , (2) 

where we normalize to as = as(2^)/(47l). The leading-order (LO) terms in Eq. Q are identical to 
the space-like case of the initial-state parton distributions, which is the so-called Gribov-Lipatov rela- 
tion |[37l . Also at the next-to-leading order (NLO), the functions Pj,- (x) are related to their space-like 
counterparts by a suitable analytic continuation, see e.g. |[38l . Following these ideas, we have derived at 
NNLO the three distinct non-singlet functions P'^^^ (with ^ = ib,v) |[39l and the diagonal singlet quan- 

titles Pqq , Pgg f40l from the corresponding space-like results pTT] and f^^l using two independent 
methods. One approach relied on mass factorization and the structure of the infrared singularities, the 
other implemented the concept of universal (kinematics independent) splitting functions as conjectured 

_ T f2) T 

in |[43l . The off-diagonal three-loop time-hke quantities Pgq , Pqg are presently still unknown. 
Theory of (extended) dihadron fragmentation functions 

Marco Radici 

Dihadron Fragmentation Functions (DiFF) describe the probability that a quark hadronizes into two 
hadrons plus anything else, i.e. the process q hAhsX. They can appear in lepton-lepton, lepton-hadron 
and hadron-hadron collisions semi-inclusively producing final-state hadrons. A thorough investigation 
of their formal properties was performed in the last decade, identifying the leading [44] and sub-leading 
twist 1 45 1 terms, and recognizing the possibility of separating the diagonal and interfering contributions 
for the {hA^hs) pair being in different relative partial waves ||46J. 

At present, the most important application deals with DiFF as analyzers of the fragmenting quark 
spin and regards the extraction of the transversity distribution in the nucleon [1471 . whose knowledge 



9 



is a basic test of QCD in tlie nonperturbative domain P8l . The extraction proceeds via the measurement 
of a single-spin asymmetry, that has been recently performed by the HERMES (49] and COMPASS f5Ql 
collaborations. Information on DiFF can be extracted by the e^e^ — > (hAihAi) {hBihBi)^ process [5l|, 
wliich is currently being measured by the BELLE collaboration |[52]| . or by models ll53]l . 

Recently, DiFF have been recognized as a crucial ingredient for canceling all coUinear singularities 
in the NLO calculation of the cross section for e^e^ —^hAh^X [54|. When DiFF depend only on the 
fractional energies of the two final hadrons, there is no way to distinguish between the usual q ^ Ha hsX 
process and the mechanism qAqs^hA hsX, where the initial parton branches into two other partons 
qA, qB, each one subsequently hadronizing in one single hadron. The net outcome is that evolution equa- 
tions for DiFF contain an inhomogeneous term involving single-hadron fragmentation functions Ii54] . 

However, most of the experimental information about DiFF consists of spectra in the invariant mass 
of the hadron pair. We keep the explicit dependence upon M| inside DiFF and we define the so called 
extended DiFF (cxtDiFF). Using the jet calculus technique, we show that this (soft) scale breaks the 
degeneracy between q HaHsX and q —>^ qAqB hAhBX mechanisms, producing for the extDiFF the 
usual DGLAP evolution ll55l . Therefore, we deduce that the NLO cross section for the e^e^ — > HaHsX 
process, when differential also in Mf^ and for the hadron pair belonging to the same jet, can be described 
as the factorized convolution of extDiFF and of the same kernel used for the single-hadron fragmentation 
case. Finally, we show that it is possible to separately study the evolution of each single term in the 
expansion in partial waves, and we give a preliminary numerical example. 

Fragmentation functions in e+e collisions at low 

David Kettler 

We report measurements of transverse momentum pj spectra for ten event multiplicity classes of 
p-p collisions at y/s = 200 GeV at STAR. By analyzing the multiplicity dependence we find that the 
spectrum shape can be decomposed into a part with amplitude proportional to multiplicity and described 
by a Levy distribution on transverse mass mj, and a part with amplitude proportional to multiplicity 
squared and described by a Gaussian distribution on transverse rapidity yj. The functional forms of 
the two parts are nearly independent of event multiplicity. The two parts can be identified with the soft 
and hard components of a two-component model of p-p collisions. This analysis then provides the first 
isolation of the hard component of the pj spectrum as a distribution of simple form on jt- 

We analyze the energy scale dependence of fragmentation functions from e^e^ collisions using con- 
ventional momentum measures Xp and and rapidity y. We find that replotting fragmentation functions 
on a normalized rapidity variable results in a compact form precisely represented by the beta distribution, 
its two parameters varying slowly and simply with parton energy scale Q. Dijet multiplicities are used 
to constrain these parameters at lower energies. The resulting parameterization enables extrapolation of 
fragmentation functions to low Q in order to describe fragment distributions at low transverse momentum 
Pt- These results compared favorably to conventional representations of fragmentation functions. 

Finally, we analyze the minimum bias two-particle correlations of p-p collisions. Cuts are used to 
isolate the jet-like peak in the angular correlations and the corresponding two-particle yr correlations are 
compared to the fragmentation function parametrization derived from e+e^ collisions. 
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Jet fragmentation at MLLA and beyond 

Redamy Perez-Ramos 

The hadronic -spectrum inside a high energy jet is determined including corrections of relative 
magnitude o(y^) with respect to the Modified Leading Logarithmic Approximation (MLLA), in the 
limiting spectrum approximation (assuming an infrared cut-off Qq = Aqcd) and beyond {Qq ^ Agco)- 
The results in the limiting spectrum approximation are found to be in impressive agreement with pre- 
liminary measurements by the CDF collaboration, unhke what occurs at MLLA, pointing out small 
overall non-perturbative contributions. Within the same framework, 2-particle correlations inside a jet 
are also predicted at Next-to-MLLA and compared to previous MLLA calculations. MLLA corrections, 
of relative magnitude o{^/U^) with respect to the leading double logarithmic approximation (DLA), 
were shown to be quite substantial for single-inclusive distributions and 2-particle correlations ||561l57i 
Therefore, it appears legitimate to wonder whether corrections of order o{^/o^s), that is next- to-next- to- 
leading or next-to-MLLA (NMLLA), are negligible or not. 

The starting point of this analysis is the MLLA evolution equation for the generating functional of 
QCD jets [58 1 . Together with the initial condition at threshold, it determines jet properties at every ener- 
gies. At high energies one can represent the solution as an expansion in Then, the leading (DLA) 
and next-to-leading (MLLA) approximations are complete. The next terms (NMLLA) are not complete 
but they include an important contribution which takes into account energy conservation and an improved 
behavior near threshold. Some results for such NMLLA terms have been studied previously for global 
observables and have been found to better account for recoil effects. They were shown to drastically af- 
fect multiplicities and particle correlations in jets: this is in particular the case in [59], which deals with 
multiplicity correlators of order 2, and in ll60l . where multiplicity correlators involving a higher number 
of partons are studied; in particular, the higher this number, the larger turn out to be NMLLA corrections. 

The present study makes use of this evolution equation to estimate NMLLA contributions to our dif- 
ferential observables. It presents the complete calculations of the single inclusive kj^ distribution leading 
to the results published in [61], and extends them to 2-particle correlations inside a high energy jet. The 
obtained agreement between the CDF results and the NMLLA distributions over the whole -range is 
excellent. The NMLLA calculation is in particular able to capture the shape of CDF spectra at every jet 
hardness Q. Conversely, predictions at MLLA prove only reliable at not too large kj^. The domain of 
validity of the predictions has been enlarged to larger k\^ (and thus to larger x since Y is fixed) computing 
from MLLA to NMLLA accuracy. This agreement further supports the Local Hadron Parton Duality 
(LPHD) hypotheses l.62il . 

The MLLA, NMLLA equations for correlators were analyzed and solved iteratively This al- 
lowed us to generalize the result previously obtained by Fong and Webber in [63] that was valid in 
the vicinity of the maximum of the single inclusive parton energy distribution ("hump"). In particu- 
lar, we have analyzed the regions of moderately small x above which the correlation becomes "nega- 
tive" {C — I < 0). This happens when suppression because of the limitation of the phase space takes 
over the positive correlation due to gluon cascading. This region turned out to be larger in NMLLA 
than in MLLA. Also, the correlation vanishes {C — > 1) when one of the partons becomes very soft 
(^ = In 1/x — > y = InEQ/Qo). The reason for that is dynamical rather than kinematical: radiation of a 
soft gluon occurs at large angles which makes the radiation coherent and thus insensitive to the internal 
parton structure of the jet ensemble. Qualitatively, our MLLA and NMLLA result agrees better with 
available OPAL data than the Fong-Webber prediction. There remains however a significant discrep- 
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ancy, markedly at very small x. In this region non-perturbative effects are likely to be more pronounced. 
They may undermine the applicability to particle correlations of the local parton-hadron duality consid- 
erations that were successful in translating parton level predictions to hadronic observations in the case 
of more inclusive single particle energy spectra ll58l . 

3 Experimental aspects of parton fragmentation in the vacuum 

Colour coherence and a comparison of the fragmentation of gluon and quark jets 

Klaus Hamacher 

Coherence is an important property of gluon radiation in the hadronic final state of inelastic inter- 
actions and a direct consequence of the quantum nature of strong interactions. Theoretical calculations 
observing coherence predict the evolution scales of jets inside of multi-hadronic events. Typically these 
scales are transverse-momentum-like and fulfill Lorentz-invariance in contrary to the often used jet en- 
ergies. A clear demonstration of colour coherence is obtained from the measurement of the ratio of the 
soft hadron multiplicity in narrow cones oriented perpendicular to the plane of a three-jet event and the 
axis of a two jet event II641I65II . This ratio is predicted as the product of the colour factor ratio Ca/Cf and 
an event scale rj^. The rj^ -dependence allows to compare the amount of coherent soft radiation due to 
the gluon and the quark colour charge. The prediction without any free parameter describes the data very 
well. A destructive l/N^ suppressed term included in rj^ is required to describe the data. The colour 
factor ratio when fitted to the data is Ca/Cf = 2.211 ib0.014(stat.) it 0.053(syst.). The reason for this 
good agreement of the leading order prediction with the data lies in the soft, large angle hadrons used as 
testing probes. These hadrons due to their large wavelength are insensitive to higher order structures and 
also to finite energy and leading particle effects. 

A precision study of the overall multiplicity of three-jet events Il66ll65]| also relies on properly de- 
fined scales in order to describe the topology dependence of the multiplicity. Two different choices are 
available as the division of the event multiplicity in a qq and a gluon part is ambiguous. The most natural 
choice describes the data whereas the choice minimising the amount of hadrons assigned to the gluon 
jet (chosen by OPAL) fails to agree. The colour factor ratio was fit to the slope of the scale dependence 
of the multiphcity yielding Ca/Cf = 2.261 ± O.OUstat. ± 0.036exp. ± 0.052theo. ± 0.041c!us.- This resuU 
together with a measurement of the P-function of the thrust [.67il which agrees to the one of the coupling 
in NLO restrict the gauge group of strong interactions to SU(3). From the measurement the multiplicity 
of gluon-gluon colour singlet events can be extracted. This multiplicity rises about twice as fast with the 
energy scale as in the qq-case clearly demonstrating the higher colour charge of the gluon. 

The precise agreement of the above discussed measurements to the colour factor expectation leaves 
little room for further differences in the particle production from gluons and quarks. The measurements 
mainly admit differences for leading particles. Such are indeed observed for baryons f68l, where an 
excess of about a factor two (i.e. ~ 2 • Ca/Cf) is observed for leading particles. This excess is under- 
stood within the string fragmentation model due to the higher number of possible splitting processes into 
baryons. Cluster models here fail if no fundamental splittings of gluons into diquark-anti-diquark pairs 
are foreseen. 

It has been found in Monte Carlo studies that gluon fragmentation functions obtained from three- 
jet events and hypothetical gluon-gluon events disagree in specific parts of phase-space [M. Tasevsky]. 



12 



Therefore, the resuhs obtained from three-jet events were quantified as "biased". Tliis distinction seems 
unfounded as tiiree-jet events form the really accessible source of data. Moreover, if the observed differ- 
ences were essential they need to be brought into accord with the QCD factorisation theorem. It turns out, 
however, that the discrepancies are consequence of hadronisation and the necessities of the measurement 
process ||69]| . The reconstruction of the parton kinematics due to confinement must proceed by a direct 
identification of the jet with the parton kinematics. However, the hadronisation transition then leads 
to a smearing of the parton kinematics. Typically this smearing vanishes oc l/Ejet- The parton energy 
inferred from the jets enters in the definition of the scaled momentum x^, = Phadmn/ Pg/q and in turn espe- 
cially strongly varying fragmentation functions are influenced by error propagation. In fact the effect is 
stronger for the gluon compared to the quark fragmentation function due to the stronger fall-off, as well 
as due to the somewhat increased hadronisation smearing for gluon jets. A numerical calculation shows 
ratios of up to a factor two between the smeared, measurable fragmentation function and the underlying 
partonic one for gluons at high Xh and small Ejgf The difference detected for gluon fragmentation [M. 
Tasevsky] can be fully explained by the hadronisation/measurement process. It is important to note that 
similar problems unavoidably exist in any measurement involving jets, e.g., when gluon fragmentation 
is measured in pp 2 jets. The hadronisation smearing cannot be properly included in an experimental 
unfolding as it will depend on the fragmentation function to be measured. Still, it is possible to use such 
measurements if the hadronisation smearing is included in the phenomenological fitting procedure of the 
fragmentation functions. 

ALEPH results on quark and gluon fragmentation 

Gerald Rudolph, for the ALEPH collaboration 

Parton fragmentation is studied in e^e^ annihilation into hadrons in the energy range 91 - 207 GeV 
using the ALEPH detector at LEP. Scaling violation in the flavour-inclusive x-distribution of charged 
particles has been observed which is, for x > 0.1, qualitatively reproduced by global NLO parametrisa- 
tions IITOI TTI]. The data are better described by the QCD Monte Carlo (MC) models. The energy variation 
of < Nch > is well explained by the QCD-MC models and by a two-parameter 3NL0 calculation. 

Inclusive x distributions of a variety of identified light mesons and baryons have been measured at 
the Z pole 117211731 . These data supplemented with pj and event shape distributions were used to tune 
the free parameters of the QCD-MC programs JETSET, ARIADNE and HERWIG [71.24]. The long- 
standing discrepancy at /Jx.out > 1 GeV/c can be slightly improved by using the new ;77-ordered parton 
shower 1751 . The fragmentation functions of D*+ and B mesons as well as the relative production rates 
of the 0^, 1^ and higher meson spin states in the c and b sector have been measured f76l 17711781 . To 
simulate Bose-Einstein correlations in hadronic Z decays, the parameters X and a of the ^£'32 model of 
PYTHIA have been determined |[79l . 

The gluon and quark fragmentation functions have been measured from symmetric 3-jet events at 
two scales using Z^-tagging techniques II80II81II . As expected, the gluon fragmentation function is steeper 
in X than the quark fragmentation function. The failure of the QCD-MC models to describe the x > 0.4 
region of the gluon fragmentation function, first noted in [18711821 . is confirmed with high statistics in 
3-jet events of general topologies with high statistics (preliminary). The rate of neutral gluon jets with 
a central rapidity gap is found to be higher than expected from JETSET or ARIADNE [,74. 1 . The excess 
is located at low invariant mass (0.8-2.2 GeV/c^, preliminary), as first observed by DELPHI ll83l . The 
colour-reconnected versions of these generators predict a much higher rate and thus can be excluded fJ4\ . 
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Colour flux studies in quark and gluon fragmentation in the DELPHI experiment at LEP 

Brigitte Buschbeck, for the DELPHI collaboration 

The occurrence of colour reconnection between two quark-antiquark (colour triplet) strings |843, 
Bose-Einstein correlations of their decay products fSF], and the possible occurrence of colour octet 
fragmentation with the formation of gluonic systems in the fragmentation of gluons f86l is studied in 
three contributions. 

The first contribution investigates colour reconnection directly in e^e~ ^ W W ^ {q\ - q2) {q3- q4), 
i.e. between the colour strings spanned by each W. The second contribution investigates Bose-Einstein 
correlations (BEC) between the decay products of the two different W bosons. Both effects are de- 
termined by the space-time overlap of the two colour flux tubes or their decay products. The data are 
compatible with some reconnection effects, however with a large statistical error. Bose-Einstein corre- 
lations between the decay products of different W's are seen by DELPHI (significance 2.4a ), however 
they show up weaker than predicted by a Monte Carlo which includes full inter W correlations. Since 
all other LEP2 experiments do not show any signal, it is argued, that BEC may be diminished or even 
absent between the hadronisation products of different strings. This can be a valuable tool to study the 
string structure in hadron-hadron and heavy ion reactions. 

The third contribution investigates the leading hadrons in gluon fragmentation in 3-jet (e^e^ qqg) 
events f86l. It shows that gluon jets produce more neutral leading systems than predicted by the LUND 
string models (JETSET, ARIADNE). This is however predicted, if gluons {g) fragment in some cases as 
a colour octet field which is then neutralised by a gg pair The effect is enhanced by demanding a rapidity 
gap and amounts at a rapidity gap = 1.5 of about 10 % . The concentration of the excess of the neutral 
systems at low invariant mass could be an indication for a hitherto undetected fragmentation mode of the 
gluon via the formation of a gluonic system. 

Fragmentation functions of quark and gluon jets as measured by OPAL 

Marek Tasevsky, for the OPAL collaboration 

Scaling violations of quark and gluon jet fragmentation functions are studied in e+e^ annihilations 
at = 91.2 and 183-209 GeV using data collected with the OPAL detector at LEP. The scale depen- 
dence of the flavour inclusive, udsc and b fragmentation functions from unbiased jets is measured at 
y/s/2 = 45.6 and 91.5-104.5 GeV. Biased jets are used to extract the flavour inclusive, udsc and b, and 
gluon fragmentation functions in the ranges Q-^^t = 4-42, 4-105 and 4-70 GeV, respectively, where Q-^et 
is the jet energy scale. Three methods are used to extract the fragmentation functions, namely the b-tug 
and energy-ordering methods for biased jets, and the hemisphere method for unbiased jets. The results 
obtained using these methods are found to be consistent with each other The udsc jet results above 
the scale of 45.6 GeV, the gluon jet results above 30 GeV (except for the scale of 40.1 GeV), and the 
b]&t results at all scales except 45.6 GeV represent new measurements. The results of this analysis are 
compared with existing lower energy e+e^ data and with previous results from DELPHI and OPAL. The 
overall consistency of the biased jet results with the unbiased jet results suggests that Qj^, is a generally 
appropriate scale in events with a general three-jet topology. The scaling violation is observed to be 
positive for lower xe and negative for higher xe, for all the types of fragmentation functions. The gluon 
jet fragmentation function exhibits stronger scaling violation than that of udsc jets. 

The bias of the procedure used to construct biased jet fragmentation functions is estimated by study- 
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ing hadron level Monte Carlo generator events. In explaining the observed differences between biased 
and unbiased jet results, we note the effects of non-negligible masses of hadrons and ft-quarks at low 
scales. Due to the considerable bias found for the gluon jet fragmentation functions in the region of 
xe > 0.6, precautions should be taken when comparing the biased gluon jet results with theory. The 
data are compared to the predictions of NLO calculations. In a wide range of scaled momentum xe, all 
calculations satisfactorily describe the data for the udsc jet fragmentation functions. The description is 
worse and the spread between the predictions larger for the b and gluon jet fragmentation functions, in 
particular in regions of very low and high xe- The data are also compared with predictions of three Monte 
Carlo models, PYTHIA 6.125, HERWIG 6.2 and ARIADNE 4.08. A reasonable agreement with data is 
observed for all models, except for high xe region with small scales (< 14 GeV) in case of the udsc and 
gluon jet fragmentation functions. The charged particle multiplicities of udsc, b and inclusive hadronic 
events are obtained by integrating the measured fragmentation functions. All values are found to be in 
agreement with previous measurements, where available. 

The first experimental study to use the jet boost algorithm, a method based on the QCD dipole model 
to extract properties of unbiased gluon jets from e+e^ —>^<iQg events, has also been presented. We test 
the jet boost algorithm using the Herwig Monte Carlo QCD simulation program, comparing the results 
of this method to those derived from unbiased gluon jets defined by hemispheres of inclusive gg events 
from a color singlet point source. We find that the results of the jet boost algorithm for the multiplicity 
distribution are in close correspondence to those of the gg hemispheres for jet energies E* larger than 
about 5 GeV. For the fragmentation functions, the results of the two methods agree to good precision for 
E* > 14 GeV. Therefore, the fragmentation function of unbiased gluon jets have been measured only at 
two points, namely at 14.24 and 17.72 GeV. 

Parton fragmentation studies in BaBar 

Fabio Anulli, for the BaBj^ collaboration 

The BAB^ detector [87], operating at the asymmetric S-factory PEP-II, has been optimized for CP 
violation studies in B meson decays. However, the high luminosity and the high detector performances, 
allow also precise measurements of different aspects of strong interactions. We present preliminary mea- 
surements of inclusive momentum spectra of a variety of light mesons and baryons, namely 7i*, K^, 
r\, and p/ p, at an energy of ^/s = 10.54 GeV, below the BB production threshold. The data cover the 
full scaled momentum range with high precision (few percent relative), allowing sensitive tests of QCD 
calculations and fragmentation models. Comparison with results from higher energies shows significant 
scaling violation at high scaled momentum values, for pions and kaons, while data for p/p do not present 
any clear difference between 10 and 90 GeV. The same spectra have been measured also in hadronic 
decays of the T(45), providing significant inputs to model the inclusive properties of B meson decays. 

We present also a measurement of the inclusive spectra of the lightest charmed baryon, the A+ f88T|. 
The scaled momentum distribution is in agreement with the existing data and shows a maximum at 
Xp PS 0.6, followed by a sharp decrease. The total measured production rate is Na^ = 0.057 it 0.002(exp) 
±0.015(BF) A+ per hadronic event, where the dominant uncertainty is given by the branching fraction 
into the reconstructed mode, A+ —>■ pK^n^. The production of A+A^ pair has been measured at a rate 
about four times higher than expected from the measured single A+ production rate; the event topology 
is not consistent with four-baryons production, while an average of about four popcorn mesons are pro- 
duced in association with the A+A^ pair. This result is compatible with long distance baryon number 
conservation, as already pointed out by the CLEO experiment 11891 . 
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The very high luminosity at the B-factory allows the study of e^e^ annihilation into exclusive, low- 
multipUcity final states. Such processes dominate at low but their cross sections fall much faster 
than ~ 1/s, and are several order of magnitudes smaller at ^/s ~ 10 GeV. Their measurement in terms 
of total cross section and internal structures, provide a rich testing ground for QCD. As an example, we 
show here the first observation of e^e^ annihilation in two hadrons via two-virtual-photon-annihilation 
(TVPA), in the processes e^e^ — > p^p" 7i+7l^7i+7i^ and e+e^ (|)p° K^K^n^n^ , whose final 
states have positive charge conjugation, C = +1, and cannot be produced by a single y [90]. Annihilation 
of e~^e^ pairs at low is studied in BaBar via initial state radiation, as shown in detail in [S. Pacetti]. 

Light-quark and spin-dependent fragmentation functions at BELLE 

RalfSeidl, for the BELLE collaboration 

The BELLE experiment fOP] at the asymmetric e^e~ collider KEKB f92l has taken a large amount 
of data which can be used for the study of fragmentation functions. In particular the more than 60 fb^ ^ 
of continuum data contain only udsc quark-antiquark pairs and no B meson decays. In leading order, 
the normalized yield of inclusive identified hadrons can be directly related to the sum of fragmentation 
functions of all accessible quarks and antiquarks. With the statistics available systematic uncertainties 
will be dominating. While the momentum scale and smearing are very well understood the amount of 
misidentified hadrons has to be corrected. For this purpose particles where the type of the final state 
particles is known from the decay D*,A are used to obtain the particle detection efficiencies and fake 
rates. Results on the unpolarized fragmentation functions are expected soon. 

A second type of fragmentation functions studied in BELLE treats the fragmentation of transversely 
polarized quarks. The most prominent of these fragmentation functions is the Collins function 1931 . 
Nonzero results from a 29 fb^' continuum data sample have been published [94 1. The statistics of these 
results have been increased by nearly a factor of 20 as for these measurements also the data taken under 
the T(45) could be included due to a vanishing B-meson background in this analysis. 

A global analysis of the published BELLE Collins asymmetries and the semi-inclusive DIS results 
from HERMES r95l and COMPASS [96] containing the Collins effect in conjunction with the quark 
trans versify distribution were successfully performed by Anselmino and collaborators [97|. The obtained 
transversity distribution seems to be smaller than expected from bounds and the Lattice QCD predictions. 
However many aspects, such as the evolution of the Collins function are still open questions. Another 
spin dependent fragmentation function, the interference fragmentation function is also being studied at 
BELLE and results are expected soon. 

Fragmentation studies at CLEO 

David Besson, for the CLEO collaboration 

Over the last two decades, CLEO has made a number of measurements testing fragmentation of 
quarks in qq events vs. digluon fragmentation vs. three-gluon fragmentation at center of mass energies 
around 10 GeV. Not reported in my talk are older results on suppression of decuplet baryons to octet 
baryons in inclusive qq production on the continuum, inclusive pion, kaon and proton fractions in 3- 
gluon production vs. qq production. Reviewed results include: 

1. The older correlated Ac/Ac enhancement observed on the continuum ||89l . now being redone with 
higher statistics by BABAR, which showed the original enhanced baryon-antibaryon production 
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2. Similar A/A correlated production (9E], which showed a similar effect, even after correcting for 
c-meson cascades. As with the previous process, JETSET 7.4 underestimates this enhancement by 
about a factor of 3. 

3. Inclusive deuteron production in gluonic decays of the narrow bb T resonances |l99l, which indi- 
cates a 3-fold enhancement of per-event deuteron production in ggg events compared to qq events. 

4. A comparison of particle production quark vs. gluon fragmentation, using qq (y) ISR events on 
the continuum to tag qq events of a particular vs. ggy events from the narrow resonances to tag 
gg fragmentation at the same Q^. After confirming the large excess of baryons in ggg production 
relative to qq production, (about double the JETSET 7.4 prediction), the photon-tagged analysis 
indicates only a modest (and less momentum-dependent) baryon enhancement in gg events, and 
better consistency with JETSET 7.4. 

5. A comparison of direct photon production in quarkonium decays. This, in turn, depends on the 
color-octet and color-singlet characteristics of the recoil gg states. Newer data on the direct photon 
momentum spectrum from the \|/ (presented first at this conference) are consistent, once accounting 
for two-body radiative decays, with older data from the T IIIOOI . This is somewhat surprising, given 
the larger relativistic corrections in the charmonium system. We also note that the rate of ggy/ggg 
is about a factor of 50% larger than the value that extrapolates directly from the T, taking into 
account the running of from the T to the \|/ mass region. 

6. A comparison of the average charged multiplicity in two-gluon vs. qq events IIIOII . which showed 
that, at our center of mass energies, the average multiplicities are consistent with unity, to within 4%. 

7. Inclusive production of charm from the Xi, states (/ = 1; decaying into qq, with a nearly on-shell 
gluon). The level of production observed (approximately 10%) is consistent with a recent 
color-octet-i-color-singlet calculation of Braaten et al. [|l02l . 

Particle production and fragmentation studies at HI 

Daniel Traynor, for the HI Collaboration 

Results from HI were presented for the average charged track multiplicity and for the normalised 
distribution of the scaled momentum, Xp, of charged final state hadrons as measured in deep-inelastic e-p 
scattering at high in the Breit frame of reference II103II . The analysis covers the range of photon vir- 
tuality 100 <Q^ < 20000 GeV^. The results are compared with e^e^ annihilation data and with various 
calculations based on perturbative QCD using different models of the hadronisation process. 

The results broadly support the concept of quark fragmentation universaUty in e-p collisions and 
e^e^ annihilation. A small multiplicity depletion compared to e^e^ is observed at low Q which can be 
attributed to higher order QCD processes occurring as part of the hard interaction in e-p scattering but 
not in e^e^ annihilation. At high Q a large depletion is also observed. In the low and high Q regions, 
where the comparison to e^e^ is poor, the Monte Carlo models are able to provide a better description 
of the data. The results are compared with NLO QCD calculations as implemented in the CYCLOPS 
program. All three parametrisations of the fragmentation functions used in this program fail to describe 
the scaling violations seen in the data. 

Results from a study of mini jets in deep-inelastic electron proton scattering were presented with 
the aim of finding evidence for hadronic activities in excess to those expected from the primary interac- 
tion II104II . The analysis is performed separately for the inclusive jet sample, which defines the towards 
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direction, and for a dijet sample where the second jet is required to have an azimuthal angle larger than 
140 degrees with respect to the leading jet, which defines the away region. The dijet sample is split into 
two samples which are enhanced in direct photon and resolved photon processes. The transverse region 
between the toward and away regions should be sensitive to additional hadronic activity. The results are 
compared to various QCD based models. The analysis covers the range 5 < < 100 GeV^, the leading 
jet is required to have a prjet > 5 GeV/c as are both dijets, the mini jets must have a. px > 3 GeV/c. 

An overall good description of the data in the "toward" and "away" regions is given by all models 
for both the inclusive sample and the dijet sample, which proves that the models are able to describe the 
primary process. In the transverse region the predictions of the models without multiple interactions are 
generally undershooting the data for the inclusive sample and for the dijet sample with resolved photon 
events. The inclusion of multiple interactions improve the agreement significantly but not completely. 

Particle production at ZEUS 

David H. Saxon, for the ZEUS collaboration 

Results are presented and compared to theory and models on a range of topics using part or all of the 
full HERA luminosity of 0.5 pb^^ ifTMl : 

• Charged Multiplicities and scaled momenta. The use of the Breit frame for comparison to e^e~ 
results and the choice of scale are explained. The benefits of using the hadronic centre of mass 
frame are shown. Scaled momentum distributions are shown for 10 < 2^ < 4098O GeV^ and 
compared to MC and MLLA predictions based in LEP fits. Scaling violation is presented. 

• Strange Particle Production. ^*',A and A distributions are presented. Compared to expectation 
there is an excess of A's in resolved photoproduction. 

• Anti-deuteron and antiproton production. Signals identified by dE /dx are shown, 'p and p rates 
are equal whereas d production is suppressed. 

• Charm fragmentation anJFj" . D*^ ,D'^ ,Ds and A+ signals are presented. Fragmentation 
parameters are extracted and are mostly consistent with e^e^ results, with an excess of A^. 

• Excited charm and charm-strange mesons. Signals are seen for a range of mesonic states. Angular 
distributions indicate the likely spins. 

• Baryons decaying to strange particles. Signals areseenfor Af.(2286),A(1520),S(1320),r*(1385) 
and other states. 

• K^Kf Bose-Einstein correlations. The formalism is recalled and comparison made to K^K^. The 
confusion caused by /o (980) production is described. 

Jets and MLLA studies at tlie Tevatron 

Alexandre Pronko, for the CDF collaboration 

The evolution of jets is driven by the emission of gluons at very small transverse momenta with re- 
spect to the jet axis. Analytical predictions are based on Next-to-Leading Log Approximation (NLLA) |[58i 
calculations supplemented with the hypothesis of Local Parton-Hadron Duality (LPHD) 162). NLLA 
provides an analytical description of parton shower formation, while LPHD states that the hadronization 
process takes place locally, and, therefore, properties of hadrons and partons at the end of the parton 
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shower are closely related. Detailed studies of jet fragmentation allow one to better understand the rel- 
ative roles of perturbative parton showering and non-perturbative hadronization in shaping the main jet 
characteristics. The Tevatron data presents a unique opportunity to verify the validity and consistency 
of the NLLA+LPHD approach on a broader range of jet energies than that available at other machines. 
Overlap of the energy scales between Tevatron and e+e^ experiments allows for a direct comparison of 
experimental results obtained in different environments. 

The inclusive momentum distributions of charged particles, j^^, in jets were measured for dijet 

events in a wide range of dijet masses, 80-600 GeV/c^ |[T06l . The analysis was done for particles in 
restricted cones around the jet direction (6c = 0.28, 0.36, 0.47). The data were found to agree with theo- 
retical calculations in the NLLA-i-LPHD framework, in the region where NLLA is indeed applicable. A 
fit of the shape of the distributions yields the NLLA cutoff scale Q^ff = 230±30 MeV and the ratio of 
charged hadrons to the number of partons produced in a jet K^ip^^q'^ = 0.56ib0.10. 

The charged particle multiplicities in quark and gluon jets, Nq and Ng, were obtained by comparing 
charged particle multiplicities in two data samples: dijet data and photon-f-jet data [107]. These two sam- 
ples have a different quark/gluon jet content (^^60% for dijet and ~20% for y+jet events), which allows 
for a measurement of the inclusive properties of gluon and quark jets. The evolution of measured charged 
particle multiplicities in gluon and quark jets with jet hardness scaling variable Q=2Ejettan{'dc /'^) agree 
well with trends predicted by recent 3NLLA calculations |108| if Qeff is set to 230 MeV |106| and 
normalization constant is let to be only free parameter. We also measured the ratio of charged particle 
multiplicities in gluon and quark jets, r = Ng/Nq. At Q = 19 GeV, we obtain r = Ng/Nq = 1.64ib0.17, in 
agreement with 3NLLA calculations II108II . 

The measurement of the two-particle momentum correlations in jets is based on dijet events with 
the invariant mass, Mjj, in the range from 66 to 563 GeV/c^ II109I . The correlation function was intro- 
duced |[63l in terms of the parameter ^ = ln{Ejet / Pimdron) and was defined as a ratio of 2- and 1 -particle 

inclusive momentum distributions: C(A^i,A^2) = D^l^)D{i,2) = + (A^i +A^2) +C2(A^i —^^lY- The 
results were obtained for charged particles within a restricted cone with an opening angle of 6^ =0.5 ra- 
dians around the jet axis. Overall, the data and theory show the same trends over the entire range of dijet 
energies. The NLLA fits 1631 to the evolution of correlation parameters c\ and C2 with jet energy scale 
allowed to extract the value of the parton shower cutoff scale Qgff = 137^Q gg MeV. The comparison to 
Monte-Carlo revealed that both Pythia Tune A and Herwig 6.5 reproduce correlation in data fairly well. 

The measurement of the intrinsic kj (transverse momenta of particles with respect to jet axis) distri- 
butions is of particular interest because it allows to probe softer particle spectra than previously measured 
observables. The analysis was based on the same dijet events that were used in the described above mea- 
surement of the two-particle momentum correlations in jets. To compare shapes of the kj distributions 
in data with the recent theoretical calculations ||56]| . the predictions were normalized to agree with data 
in the bin —0.2 < Inikj) < 0.0. Within the validity region f56l, the agreement between data and theory 
is fairly good. The inclusion of the higher order corrections in theoretical calculations [61] improves the 
agreement over the wider range of ln{kT)- 
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4 Strange and heavy quarks fragmentation: theory and experiment 



Evidence for power corrections in heavy quark fragmentation in e^e collisions 

Matteo Cacciari 

I report on the findings of IIIIOII . whiicii compare QCD tlieoretical predictions for heavy flavoured 
meson fragmentation spectra in e^e^ annihilation with data from CLEO Ullll . BELLE HI 121 and 
ALEPH |76|. Several ingredients are included in the calculation: next-to-leading order initial conditions, 
evolution and coefficient functions; soft-gluon resummation at next-to-leading-log accuracy; a matching 
condition for the crossing of the bottom threshold in evolution II113I implemented at next-to-leading or- 
der accuracy; and important initial-state electromagnetic radiation effects in the CLEO and BELLE data. 

It is found that with reasonably simple choices of a non-perturbative correction to the fixed-order 
initial condition for the evolution, the data from CLEO and BELLE can be fitted with remarkable accu- 
racy. The main result of the analysis is however that the fitted fragmentation function, when evolved to 
LEP energies, does not reproduce fully well the D* fragmentation spectrum measured by ALEPH. Large 
non-perturbative power corrections to the coefficient functions of the meson spectrum are speculated to 
be needed in order to reconcile CLEO/BELLE and ALEPH results. The data do not allow, however, 
to distinguish between a standard C/Q^ correction with a fairly large coefficient (C ~ 5 GeV^) or an 
unexpected linear term, C/Q, with a more reasonable coefficient (C ~ 0.5 GeV). 

Heavy quark fragmentation with an effective coupling constant 

Gennaro Corcella 

The transition of partons into hadrons, for the time being, cannot be calculated from first principles 
QCD, but is usually described in terms of phenomenological models, containing parameters which are 
to be fitted to experimental data. Alternatively, one can model non-perturbative corrections by means of 
an effective strong coupling constant, which does not exhibit the Landau pole any longer and includes 
absorptive effects due to parton branching lill4J . 

The talk discusses bottom- and charm-quark fragmentation in e^e^ annihilation, using the effective- 
coupling model to include hadronization effects. Perturbative quark production is treated in the frame- 
work of perturbative fragmentation functions [115], which factorizes the production of a heavy quark 
as the convolution of a coefficient function, describing the emission off a massless parton, and a per- 
turbative fragmentation function, associated with the transition of a light parton into a massive quark. 
The perturbative fragmentation function follows the DGLAP evolution equations, whose solution allows 
one to resum the large logarithm of the heavy-quark mass exhibited by the fixed-order calculation. In 
our computation, we use NLO MS coefficient functions and NLO initial condition of the perturbative 
fragmentation function. The DGLAP evolution equations are implemented in the non-singlet sector, to 
NLL accuracy, and threshold resummation is performed in the NNLL approximation, in both coefficient 
function and initial condition. The effective coupling constant, evaluated to NNLO accuracy, is the only 
source of non-perturbative corrections. In order to account for the theoretical uncertainty on our pre- 
dictions, we vary the parameters entering in the perturbative computation, such as renormalization and 
factorization scales, the ^-quark mass and as{niz), within conventional ranges. 

Following 111161 . we consider data on S-hadron energy distributions at the Z" pole, collected by 
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LEP and SLD experiments. In particular, the ALEPH collaboration reconstructed only B mesons, while 
OPAL and SLD had also a small fraction of Zj-flavoured hadrons. We find that our calculation, provided 
with the effective-coupling model, is able to give a reasonable description of the considered data, within 
the experimental and theoretical uncertainties. We also compare our model with data in Mellin moment 
space, measured by the DELPHI collaboration, and obtain fair agreement with the first five experimental 
moments, within the errors. We then investigate data on charmed hadrons from ALEPH and B-factory 
experiments, along the lines of [ 1 17 1. Our result is that the effective-coupling model is able to acceptably 
reproduce the ALEPH D*^ spectrum, while serious discrepancies are present when comparing with D 
and D* data from CLEO and BELLE. In fact, other analyses had previously shown that, even with a 
flexible non-perturbative fragmentation function, provided with three tunable parameters, one was not 
able to reproduce at the same time both ALEPH and B-factory D-data. Nonetheless, we still manage to 
describe quite well the first ten Mellin moments of all considered data from ALEPH, CLEO and BELLE. 

Towards a model independent approach to kaon fragmentation functions 

Ekaterina Christov^ 

We derive different measurable quantities in kaon production in e^e^ ^K+X and in Semi-inclusive 
DIS (SIDIS) eN ^ eK + X that measure certain Non singlet (NS) combinations of quark fragmentation 
functions (FF) for kaons I1II8.I . These quantities are obtained in a model independent way and hold in 
any order of QCD. In SIDIS, if charged K"^ are measured we have, of/^^^ = G§^ — a§ : 

^K^-K- ^ (4uy(S)Du + dv(E)Ddf^-'^' 0il + ia,/{2n))Cgg) (3) 
of-"' ~ (uv+dv) 0(1 + (a,/(27i))C,,) ® (4A, + 0^)^^-^" (4) 

that present two measurables for the two unknown EEs: D^^^^ and D^^^^ . In particularly, this 
allows one to test the usually made assumption D^^^^ = 0. Recently very precise data on unpolarized 
SIDIS was presented by HERMES II119II . If both and are measured, then the combination + 
K — 2K^, both in e^e^ and SIDIS kaon production, always measures the same NS combination of EEs 

{Du - Dd f^-^- , (a^^+^--2^ = a^^ + _ 2a^ , ^ = ^ + ^): 

j^K^+K--2Kl> = 6ao{el - ej){l + a,Q®)Df (5) 
daf^^'^^^ = [{4u-d) (1 +a.vQ^®) +a.vg®Q^®]Z)f!j-^" (6) 
daf = [(« + J) ( 1 + a, Cgg(g))+ a, g ® j"^" (7) 

These results can be used to determine the kaon EEs. Being model independent, they allow also one to 
test the parametrizations already obtained, but at different assumptions. In addition, eqs. ©-Q allow to 
compare the EEs obtained in e^e^ at rather high ~ m^, eq. dU, with those of SIDIS at quite low Q^, 
eqs. ^ and Q. Note that as D^^^^ is a NS, its evolution does not involve any other EEs. 

Heavy Quark Fragmentation Functions in e+e anniliilation at = 10.6 GeV 

Rolf Seuster, for the BELLE collaboration 

Over the years, QCD was found to be a reliable theory for many predictions. It is however restricted 
to regions of phase space where the scale of the hard interaction was (much) larger than any other scale 
in the particular event. In pushing the limit to its extreme, one learns more about the theory, its failures 

^Work supported by I3HP EU network. 
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and its successes. One limit that can be challenged at the very successfully operated Z?-factories is the 
measurement of the heavy quark fragmentation function. At Z7-factories, charm quarks are the heaviest 
quarks for which such measurements can be done, as Zj's are produced at threshold. For the measure- 
ment presented here, the scaled momentum Xp = \pcandidate\/\Pmax\ is used throughout. The BELLE 
detector jOTl is an experiment operated at the asymmetric e+e^ collider KEKB ||92l at a center of mass 
energy of about 10.6 GeV. For this analysis, about 45 • 10^ events below and 352 • 10^ events at the T(4S) 
resonance were used, corresponding to 15 fb^^ and 88 fb^' of integrated luminosity. As the contribu- 
tion from B hadron decay will only occur below xp = 0.5, combining on- and off -resonance data above 
this threshold decreases significantly the statistical uncertainty. It allows for a very precise and detailed 
description of the FF of charmed quarks into charmed hadrons, as shown in II112I . The difference of 
the on- and off-resonance data below xp = 0.5 is solemnly due to the contribution from B decays. The 
numbers reported in II112II are in good agreement with the world average. These high precision data 
will have to be tested against Monte Carlo generators. Such a comparison will improve the theoretical 
understanding of the underlying physics, as well as the description of the data by the MC itself. One 
example for such comparisons is the determination of how good commonly used EEs describe the data 
and what parameters of these functions to chose. Via a reweighting technique, special MC events were 
generated at many parameter points for many fragmentation functions. In general, the best agreement 
was found for the Bowler and the Lund EEs. The parameters at which these functions showed the best 
agreement can be found in [1 12]. The commonly used Peterson et al. fragmentation function showed by 
far the overall worst performance, not agreeing with the high precision data at all. Eor details, see Iill2ll . 

Another example, for which such comparisons could be made for are ratios of fragmentation func- 
tions. In the ratio, their systematic uncertainties will partially cancel. Whereas for the Xp dependent ratio 
of Xp(D*+)/xp(D+), the Monte Carlo can be tuned to give a perfect description of the data, the ratios 
Xp(Ds)/xp(D+) and Ac)/xp(D+) are completely mis-modelled in the Monte Carlo. No parameter will 
improve the description significantly. No conclusive reason for this discrepancy is known. 

Charmed particles production in e+e ^ cc at 10.6 GeV 

Simone Pacetti, for the BaBm^ collaboration 

During the last few years many new charmed mesons and baryons have been discovered. Some of 
these, especially in the charmed-strange mesons sector, were unexpected. The theoretical predictions 
based on potential models and the heavy quark effective theory (HQET) 11201 fail in describing mesons 
Uke the D*y(2317)+, D,/(2460)+, D*y(2860)+and X(2690)+. We focused on the analyses of such states: 

• D*^(2317)+ and D,7(2460)+ have been observed in the decays D*^(2317)+^ D+7l° [[ml and 
Dsj{2A60)^ ^ D^y, D+7i"y, D+ti+tc^ II121II . These analyses have been realized over the data col- 
lected by the BaB/U? experiment in inclusive reproduction at 10.6 GeV, with an integrated luminos- 
ity of 232 fb^^ Masses, widths, and decay modes of these mesons do not fit the HQET predictions 
based on the cs structure, hence there are speculations that both D*j(2317)+ and Dsj{2460)~^ are 
some type of exotic state, such as four-quark or molecular states II122L 

• D*^(2860)+ and X(2690)+ have been observed in their decay into D°^+^K+y with D° K-TI+, 

/r^Ti+Tl", and D+ /^"7i+7i+ II123I . This analysis is realized in 240 fb"^ of data collected by 
the BaB^ experiment. The evidence of the lightest structure is less significant, however it is 
compatible with a similar structure observed by BELLE lil24J . 

Before the advent of b factories the spectrum of charmed baryons comprised only 12 states II125II . In 
the last years BAB^ found 4 new charmed baryons and confirmed 3 other states observed by BELLE. 
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• The is a ess baryon with = 3/2+ , it is the first excited state of the Q.^ family. It has been ob- 
served by BaBar in the radiative decay: Q.* D-cJ, in cc production with an integrated luminosity 
of 231 fb^^ |126|. The mass difference with respect to the ground state Q.c is in agreement with 
lattice calculations II 12711 and theoretical expectations II128I1 . 

• The Ac(2940)+ is a udc baryon, it has been observed by BaBar decaying in D^p fl29]. It is the 
first observation of a charmed baryon which decays in a light baryon, i.e.: the proton. The data for 
this analysis refer to an integrated luminosity of 287 fb^^ 

• Two new charmed-strange baryons Ec.(3055)+ and Sc(3123)+, belonging to the g^c-family {q = 
u,d) have been discovered by the BaBar experiment, studying the A+/r^7l+ system, obtained in 
inclusive cc production at 10.6 GeV with 384 fb~^ II130I . The decay in the A+^^7i+ final state 
goes trough different intermediate states, namely: Sc(3055)"'" — > Lc{2455)^^ A+Ti+ZT^ and 
Sc(3123)+ rc(2520)++A'^ A+ti+A'^. No neutral partners have been observed. 

ISR Physics at BaBar 

Simone Pacetti, for the BaBar collaboration 

Using the initial state radiation (ISR) technique in fixed center-of-mass energy machines, like the 
flavour factories, we may mimic a typical e+e^ annihilation experiment with an energy scan in the 
center of mass. At the Born approximation, in the differential cross section for processes like e^e^ — > 
YiSR+hadrons, the function describing the radiation of a photon (Yisr) by one of the initial leptons and the 
hadronic cross section factorizes. This allows us to perform precise measurements of hadronic cross sec- 
tions at low energy (from threshold up to 4-5 GeV). The knowledge of /? = G{e^e^ —>■ hadrons)/a(e+e^ — 
iu^iu~ ) is crucial to compute the hadronic contributions to the muon anomalous magnetic moment H31I 
(low-energy values), and to the running electromagnetic coupling ttgED II132II (high-energy values). In 
addition we perform a better parameter determination for many p-, CO-, and (|)-excited states (J^^ = 1 ). 

• The light-meson final states measured by the BaBar experiment in the 1-3 GeV region are: ti+ti^ji", 
2(71+71-), K+K-%+iz- ini, 3(71+71-), 2(7i+7i-7i"), K+K-li-n+TZ") lfT34]l . 

• A new narrow (F ~ 60 MeV) vector meson called X(2175) has been discovered decaying in 
(^/o(980) HM and confirmed in the (^rj channel 11136111371 . 

• We have analyzed the reaction e+e- KKn II137II in two samples: K^K^n^ and ^'^^'+71+. The 
study of the KsK'^Tl^ asymmetric Dalitz plot allowed, for the first time, the identification of the 
two isospin components and their relative phase. We performed a global fit, exploiting five dif- 
ferent data sets, to describe the total KKn cross section in terms of intermediate resonances. We 
determined parameters for (|)- and p-excited states, with a special care to the OZI-suppressed <^n^. 
The (t)7i*' cross section, with a peak value of ~ 0.2 nb, shows two resonant structures: a broad bump 
at 1.6 GeV, compatible with a p recurrence [125] and a narrower resonance at 1.9 GeV, whose 
parameters are in agreement with those of the "dip" observed in multi-hadronic final states 111381 . 

• The BaBar experiment measured with unprecedented accuracy the differential cross section for 
e+g- pp [139]. In particular, by studying the angular distribution, BaBar provides also a 
time-like measurement for the ratio \G'e/^%\ °f ^^'^ electric and magnetic proton form factors. 
We observed, for the first time, a displacement from the "scaling behavior": \G^/Glf \ ~ 1, by 
measuring > \G'^\ near threshold. The angular distribution provides also a powerful tool to 
estimate the two-photon exchange contribution in the annihilation e^e^ — > pp. Using the BaBar 
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data this contribution turns out to be compatible with zero II140I . In addition, the pp cross section 
has two very interesting features: a non- vanishing threshold value, which is however expected as 
a consequence of the pp Coulomb interaction, and the presence of structures ("steps"), not yet 
understood, at ~ 2.2 and ~ 3.0 GeV. 

• Other baryon-antibaryon final states like: AA, and AlP have been measured by BAB^ 

with the ISR technique. In all these cases the cross section data have finite threshold values, as in 
the p'p cross section, which are not well understood. Indeed, having always neutral baryons, no 
Coulomb interaction effect can be used to cancel the vanishing of the phase space. 

Study of charm fragmentation in DIS at HI 

Zuzana Rurikovd, for the HI collaboration 

The process of charm fragmentation is studied in deep-inelastic D*^ meson production as measured 
by the HI detector at HERA |142|. Two different methods are used for the definition of the observ- 
able approximating the fraction of the 4-momentum of the D** meson with respect to the charm quark. 
The c-quark momentum is approximated in the 7* ;?-rest-frame either by the momentum of the jet includ- 
ing the D*^ meson or by the momentum of a suitably defined hemisphere which includes the D*^ meson. 

The visible phase space of this analysis is given by cuts on the photon virtuality 2 < 2^ < 100 GeV^, 
event inelasticity 0.05 <y<0.1 andD*"*" meson phase space 1.5 < pt:{D*^) < 15GeV/c, |ri(D*''')| < 1.5. 
In addition a D** jet with E^> 2> GeV in y*p-rest-frame is required in order to have hard scale for the 
events. The fractional energy distributions of D* mesons, corrected for detector effects, are used to ex- 
tract parameters for the Kartvelishvili and Peterson fragmentation functions within the framework of 
leading order -1- parton shower QCD models RAPGAP 3.1 and CASCADE 1.2 with string fragmentation 
and particle decays as implemented in PYTHIA, as well as for the next-to-leading order QCD calculation 
in the fixed flavour scheme HVQDIS fT43l . 

The fragmentation parameters extracted using both observables are in good agreement with each 
other. Furthermore the value of Peterson parameter £, extracted for the PYTHIA parameter setting II 14411 . 
is in agreement the ALEPH tuned value of 8 = 0.04, supporting the hypothesis of fragmentation uni- 
versality between e^e^ and e-p. Finally, the fragmentation of charm produced close to the kinematic 
threshold is studied. The description of this sample by the QCD models and NLO calculation is not 
good, and the extracted fragmentation parameters are significantly different from the parameters fitted 
to the nominal sample with a hard scale. This can be interpreted as a hint of inadequacies of the QCD 
models and possibly of insufficient flexibility of the simple parametrisations used for the fragmentation 
function in the phase space region close to kinematic threshold. 

Inclusive K'^K^ resonance production in e-p collisions at HERA 

Changyi Zhou, for the ZEUS Collaboration 

The lightest glueball is predicted to have j'^^ = 0+^ II1451I146I and a mass in the range 1450-1750 
MeV/c^. Thus, it can mix with qq states from the scalar meson nonet, which have 7 = and similar 
masses. In the literature, the state /o(1710) is frequently considered to be a state with a possible glueball 
or tetra-quark composition (see for example the reviews 1,147 ,,,148,1 ). 

Inclusive K^K^ production in e-p collisions dominated by photoproduction with exchanged photon 
virtuality, Q^, below 1 GeV^, at HERA has been studied with the ZEUS detector using an integrated 
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luminosity of 0.5 fb~ from HERA-I and HERA-II combined. mesons were identified tlirough the 
charged-decay mode, — > 7i+7i~, and reconstructed with various track quality cuts for selection. The 
K^K^ invariant mass distribution was reconstructed by combining two K'^ candidates. Enhancements in 
the mass spectrum have been observed attributed to the production of /2(1270)/tZ2(1320), /2(1525) and 
/o(1710). The three states /2(1270), a^(1320) and /2(1525) are all of 7^ = 2+ and so their interference is 
seen in the total cross-section. The intensity is the modulus-squared of the sum of these three amplitudes 
plus the incoherent addition of /o(1710) and the non-resonant background. The amplitudes for /2(1270), 
02(1320) and /2(1525) production were fixed at the SU(3) ratios of 5 : —3 : 2 as expected for production 
via an electromagnetic process II1491I150I . Very competitive measurements on peak position and width 
for/2(1525) and/o(1710) are done with interference fit and the overall fit describes the data very well. 
Complete systematic checks have been performed. The final values with statistical and systematical 
uncertainties were compared well with the PDG values II125I . 

RHIC constraints on fragmentation functions for strange hadrons in p-p collisions at^ = 200 GeV 

Mark Heinz, far the STAR collaboration 

Perturbative QCD has proven successful in describing inclusive hadron production in elementary 
collisions. Within the theory's range of applicability, calculations at next-to-leading order (NLO) have 
produced accurate calculations for the pr spectra of charged hadrons in different collision systems and 
energies, which has lead to claims of universality of the underlying fragmentation functions (FFs) [15l. 
With the new mid-rapidity p-p data at ^/s = 200 GeV collected at RHIC, comparisons to pQCD can now 
be extended to identified strange baryons and mesons |24|. In the last 5 years significant improvements 
have been made in the field of NLO FF studies. Several groups have updated their parameterizations 
to include not just e+e^ but also (SI)DIS and now p-p data, thus improving the constraints on the FF 
parameters. In particular the gluon-to-hadron FF was never well constrained by e^e^ and DIS data due 
to the low probability of gluon-jet production. Conversely, at RHIC we are probing a very gluon-jet rich 
environment and the majority of final states are produced by gluons. In particular baryons are dominated 
by gluon fragmentation. According to AKK, 90% (10%) of the protons below pj = 10 GeV/c are pro- 
duced from gluons (quarks). On the other hand, pions are produced in equal parts by gluons and quarks. 

Experimentally this fact can be exploited to measure the baryon-to-meson ratios and compare them 
to LO/NLO calculations. Results by the STAR collaboration confirm that LO models are not able to 
accurately describe this ratio 12211241 and underpredict the amount of baryons at low pj. For strange 
baryon production (A) the first NLO predictions equally fail to describe the data |[T2l . In an attempt to 
solve this shortcoming, one group (AKK) [4J readjusted the gluon FF to fit the data, finding = Dg/3. 
Another group (DSS) has for the first time performed separate fits to protons and anti -protons based on 
STAR data ifTTI . and provided new FFs for pions, kaons and charged hadrons from a global analysis ifTOl . 
Another subject which did not receive sufficient attention up until recently are error estimations on the 
FF parameterizations. One group (HKNS) [14] has now published a global (Hessian) error analysis for 
charged hadrons and finds that the 4 most ubiquitous FF agree within errors. However this analysis does 
not include p-p data. The STAR experiment has also measured identified particle spectra in terms of 

mj = \l Pt + motivated by earlier studies II151I . After arbitrary normalization of the different 
particle species, this representation clearly shows different shapes of the spectra for the baryons and 
mesons for mj > 2 GeV/c^. Interestingly, a simulation of the same observable with PYTHIA (LO 
pQCD) shows a similar feature when selecting gluon-jet enriched events 1,75.1 . Finally, efforts at RHIC 
are now being focused in the extraction of medium modified FFs. Several approaches are being pursued 
in parallel, such as di-hadron correlations, 3-particle correlations and also full jet reconstruction. 
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5 Parton fragmentation in cold nuclear matter: theory and experiment 



/'r -broadening in nuclear DIS 

Hans-Jurgen Pimer 

My work deals with modifications of fragmentation processes in cold nuclear matter and in the hot 
quark gluon plasma. In the first part a phenomenological analysis of deep-inelastic hadron production 
in nuclei is performed. The main theoretical picture is a three-stage model which includes partonic 
propagation, prehadron formation and absorption, and finally hadron formation. 

This model has been successfully applied to describe the HERMES data presented in [E. Aschenauer]. 
The multiplicity ratios of ll,K,p,p produced in e-Ne, He, Kr, Xe collisions as a function of the mo- 
mentum fraction of the hadrons Zh, and the photon energy v can be modelled by a prehadron-nucleon 
absorptive cross section of one third of the hadron cross section. No (2^-dependence of the multiplicity 
ratios is calculated as observed in the data. 

The second part emphasizes the /jj- -broadening in the hadron spectra as a tool to monitor the path 
length of the quark in the nucleus. The meson-nucleon interactions have a small elastic cross sections 
relative to the total cross section, and can be neglected for pr-broadening. Again the calculated differen- 
tial Apj broadening as function of Zh,'^,Q^ are compared with data. The observed log(Q^) dependence 
indicates the interleaving of parton radiation and parton-medium scattering in the evolution process. 

Both effects can be described with a new equation where a splitting kernel and a scattering kernel are 
present in the evolution equation. Consequences of this equation for the jet profiles in LHC experiment 
have been shown and further details will be given in [S. Domdey]. 

Hadronization of pions and kaons from nuclei using DIS 

Kenneth Hicks, for the CIAS Collaboration 

The CLAS experiment eg2 was run with a variety of nuclear targets using a 5.5 GeV electron beam 
from the continuous electron beam accelerator facility (CEBAF) at Jefferson Lab. The goal of this ex- 
periment is to measure observables related to the propagation of a quark (struck by the virtual photon 
from deep inelastic scattering) through cold nuclear matter. These results could be contrasted with quark 
propagation through hot QCD matter, as in the quark-gluon plasma expected to be formed in relativis- 
tic heavy ion collisions (RHIC). Quarks from RHIC are tagged by back-to-back "jets" of high energy 
hadrons formed in hard parton collisions. In contrast with RHIC results, where the quark propagation is 
severely damped by passage through hot QCD matter, the preliminary DIS results suggest that the struck 
quark propagates relatively freely (compared with fully formed mesons) through cold nuclear matter. 

The measured quantities of the DIS measurements are the squared four-momentum transfer (Q^) and 
energy transfer (v) from the scattered electron and the energy fraction (z = Eijv) and transverse mo- 
mentum (pt) of the leading hadron. By comparing the z and pj distributions of pions and kaons from 
a deuterium target with that from various nuclear targets, the quark propagation and formation into a 
hadron can be inferred. Theoretical predictions in a given model provide guidance for interpretation of 
the results. Ratios of the number of hadrons detected, normalized to the number of DIS events in a given 
kinematic bin of and v, for a nuclear target divided by the same quantities for a deuterium target are 
one quantitative measure of the effects of propagation through cold nuclear matter. The subtraction of 
the average value of p^ seen for a nuclear target minus that for a deuterium target is another measure 



26 



which, according to theoretical models, is sensitive to gluonic radiation by the quark before it forms into 
a hadron. By comparing these two quantitative measures with theoretical calculations over a variety of 
kinematics, then models of hadronization can be tested. The HERMES experiment has already published 
the nuclear attenuation ratios for several types of hadrons at DIS kinematics in the range 8 < V < 23 GeV. 
The CLAS data are for a lower range, with 2.6 < V < 4.3 GeV corresponding to shorter formation times. 
At HERMES kinematics, calculations suggest that the quark propagates through the full nucleus, fol- 
lowed by hadron formation outside of the nuclear radius. At CLAS, we expect to see hadron formation 
inside the radius of larger nuclei. 

Preliminary results for positive pions produced at CLAS for the eg2 experiment indicate a plateau for 
the difference ^pj as the nuclear size increases (targets C to Pb). This suggests that, at these kinematics, 
the average hadron formation length is shorter than the radius of Pb. In contrast, nuclear attenuation 
ratios cannot be explained by theoretical models unless some finite propagation of a "pre-hadron" (with 
a smaller cross section than fully-formed hadrons) is assumed. Attenuation ratios for neutral kaons 
appear to be smaller than those for 7i+ particles, but results are still preliminary. In conclusion, the 
CLAS data provide a new kinematic window for hadronization at the expected length scale near to the 
radius of heavy nuclei. Statistical uncertainties are small due to an experimental luminosity of about 100 
more than for the HERMES data. Both 7i+ and nuclear attenuation ratios and Ap^ are expected to be 
submitted for publication before the end of 2008. 

Characterising hadronization with nuclear DIS 

Valeria Muccifora 

A review on the recent progress in the study of the parton propagation, interaction and hadroniza- 
tion in cold nuclear matter has been presented. It has been pointed out that the cleanest environment 
to address nuclear modification of hadron production is the nuclear Deep Inelastic Scatteringdt allows 
to experimentally control many kinematic variables; the nuclear medium (i.e., the nucleus itself) is well 
known; the hadron multiplicity in the final state is low, allowing for precise measurements. Moreover, 
the nucleons act as femtometer-scale detectors of the hadronizing quark, allowing to study its space-time 
evolution into the observed hadron. 

The experimental highlights from deep inelastic lepton nuclear scattering have been shown, in par- 
ticular results from the HERMES experiment have been discussed. At HERMES, nuclear semi-inclusive 
deep-inelastic scattering is used to study the quark propagation and hadronization via single hadron mul- 
tiplicity ratio fl521, double hadron production 1 153 1 and broadening of the hadron transverse momentum 
Pt I[T54 I. While hadron multiplicity ratio is sensitive to all stages of the hadronization process, i.e. to 
the production time tp of the colorless qq system and to the formation time tf of the final hadron, the 
-broadening has been shown to be sensitive to tp because, once the qq is formed (t > tp) , no further 
broadening occurs as inelastic interaction is suppressed. 

The most recent theoretical frameworks for describing the interaction of energetic partons and space- 
time evolution of the hadronization process in nuclear DIS have been discussed. In particular models 
based mainly on parton energy loss effect II155II1561 have been compared with models based mainly 
on the interaction of the qq system in the medium and that are sensitive to the time evolution of the 
hadronization in the medium Ill571ll58t . It has been pointed out the importance to study parton propaga- 
tion and interaction in cold nuclear matter for the interpretation of data in A- A collisions. In particular 
the connections between the kinematic in DIS and in nucleon-nucleon collisions have been discussed 
and the equivalent phase spaces covered by different experiments have been compared II159II . 
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Parton propagation and hadron formation: present status, future prospects 



Will Brooks 



The process of hadronization is unique to QCD, reflecting the non-Abelian nature of the strong inter- 
action. For pragmatic descriptions of scattering phenomena, hadronization is captured by the empirically 
determined fragmentation functions. Until recently, study of the parton propagation and hadron forma- 
tion processes was limited to asymptotic properties of the final state, such as hadron multiplicities; any 
microscopic information on the space-time development of the process on the femtometer scale was 
lost. The microscopic space-time development of hadronization can now be accessed using the nuclear 
medium as an analyzer. The interactions with the nuclear medium can be used to refine the understanding 
of the mechanisms at work at 10~^^ m, yielding new insights into such fundamental processes as gluon 
emission by a quark and the enforcement of color confinement through hadron formation. Extraction of 
confinement quantities such as the lifetime of the deconfined quark now appears to be possible. Cross 
comparison of new data in a variety of reactions will potentially provide an exciting, coherent picture, 
although several important points concerning the interaction with the medium still need to be clarified. 
Foremost among these for the semi-inclusive DIS data is the question of whether the observed hadron 
attenuation is primarily due to prehadron absorption, or due to medium-simulated gluon emission. Fur- 
ther, the extent to which information derived from cold nuclear matter can be used to interpret data from 
hot nuclear matter continues to be a subject of discussion. In particular, medium-induced gluon emission 
is broadly believed to be involved in both cases, offering the possibiUty of constraining the interpretation 
of the high-energy heavy-ion data using the well-understood properties of cold nuclear matter. 

In addition to semi-inclusive DIS, Drell-Yan data offer another venue in which to study propagating 
quarks in cold nuclear matter, without the additional complication of hadron formation in the final state. 
In this talk, transverse momentum broadening in Drell-Yan and semi-inclusive DIS were compared, with 
significant differences between the Drell-Yan data from HERMES [E. Aschenauer,V. Muccifora] and 
Jefferson Lab [K. Hicks]. While several possibly relevant aspects of these differences are understood, a 
comprehensive picture awaits further theoretical development. It should be noted that the extremely high 
precision of both of the latter two, quite new, data sets is unprecedented and these data should stimulate 
a new wave of theoretical activity. The JLab data for positive pions comes as a fully three-dimensional 
function of Q^,v, and z, while the HERMES results include both 7t^ and an exploratory look at kaons. 
Further, the increase of Apj- in the JLab data appears to saturate at high nuclear mass. This is due to the 
lifetime of the deconfined quark being shorter than the dimensions of the heaviest nuclei, and a quantita- 
tive extraction of this lifetime within model assumptions is clearly feasible. An inter-comparison of the 
HERMES and JLab results for hadron attenuation shows a good consistency. The landmark HERMES 
data set has more than half a dozen different hadrons, extending to much higher v than the JLab data, and 
thus provides the best constraint on the v and flavor dependence. The JLab data, while limited thus far to 
pions and the first study, nonetheless offers two orders of magnitude more luminosity, thus permitting 
three dimensional binning of the multiphcity ratio for, e.g., Ji+. Studies of these data thus far do reveal 
systematic variations of the multipUcity ratio in Q^, v, p^, and z, with the largest variation in z and in p^, 
in agreement with HERMES. 

While further studies of the existing data are ongoing, the next significant advance in the field will 
come with the 12-GeV Jlab upgrade. The upgrade of CLAS will permit yet another order of magnitude 

increase in luminosity, adding baryon fragmentation studies, as well as a broader range in v, Q^, and z. 
The increased luminosity allows access to formation of rarer and heavier mesons such as the (|). In total, 
eleven mesons and eight baryons are included in the existing plans with an 1 1 GeV beam. 
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6 Parton fragmentation in hot & dense QCD matter: theory & experiment 



Medium-evolved fragmentation functions 

Carlos Salgado 

Medium-induced gluon radiation is usually identified as the dominant dynamical meclianism under- 
ling the jet quenching phenomenon observed in heavy-ion collisions (see e.g. II 16011 for a recent review on 
the calculation of the spectrum and the corresponding phenomenological implementation). In its present 
implementation, multiple medium-induced gluon emissions are assumed to be independent, leading, in 
the eikonal approximation, to a Poisson distribution [161.1 . In 1.162.1 we have introduced a medium term 
in the splitting probabilities 

P'°\z)=P'"'{z)+AP{z,t), (8) 

so that both medium and vacuum contributions are included on the same footing in a DGLAP approach. 
In eq. (Ull, AP{z,t) is taken directly from the medium-induced gluon radiation spectrum 

AP{z,t)^ — ——, (9) 
dzdt 



computed in the multiple scattering approximation [163]. The improvements include energy-momentum 
conservation at each individual splitting, medium-modified virtuality evolution and a coherent imple- 
mentation of vacuum and medium splitting probabilities. Noticeably, the usual formalism is recovered 
when the virtuality and the energy of the parton are very large. This leads to a similar description of the 
suppression observed in heavy-ion collisions with values of the transport coefficient of the same order as 
those obtained using the quenching weights. In a previous publication 1 164 1 it has been found that this 
formalism would lead to non-trivial angular dependences of the jet profiles when kinematic constrains 
are imposed. Similar structures were found at RHIC in two-particle correlations in the away side II165II . 

Hadronic composition as a characteristics of jet quenching at the LHC 

Sebastian Sapeta 

There are several mechanisms which may lead to modification of the hadronic composition of jets 
(hadrochemistry), e.g. exchange of quantum numbers (color, flavor, baryon number) between projectile 
and the medium, recombination of partons from the jet with partons from the medium as well as recoil 
effects. All the above may require serious modeling. Instead, we consider the framework which takes 
into account only the exchange of momentum between the medium and the developing partonic cascade. 

To calculate single particle distributions of identified hadrons we use the framework of Modified 
Leading Logarithmic Approximation (MLLA) II166I . This perturbative approach combined with the 
hypothesis of Local Parton-Hadron Duality (LPHD) was shown to reproduce correctly not only the dis- 
tributions of all charged particles but also the spectra of identified hadrons such as pions, kaons and 
protons II621I1671 . To model the medium-modification of jets we supplement the above formalism by the 
formulation of parton energy loss proposed in II168I . The effects of medium-induced gluon radiation are 
introduced by enhancing the singular parts of the leading order splitting functions by the factor 1 +/med- 
This leads to the softening of hadron spectra. The model accounts for the nuclear modification factor at 
RHIC when /med is of the order of 1. 

Our main result are the ratios of ^^/ti^ and /j^/ti* II169I . We observe a significant difference 
between the ratios in vacuum and medium jets. According to our analysis, in the presence of medium 
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the ratios increase up to 50% for kaons and 100% for protons. This enhancement seems to be a generic 
feature for radiative energy loss jet quenching models. In addition we observe only mild dependence 
on the the jet energy and the jet opening angle. The medium modified hadron spectra may be also 
superimposed on the high multiplicity environment of heavy ion collisions expected at the LHC. We 
have checked that due to the characteristically different hadrochemistry of the jet and the background 
our result concerning the ratios remains virtually unchanged. Hence, we conclude that the hadrochemical 
composition of jets may be very fragile to the medium effects and could be used as an additional handle 
to study microscopic mechanisms underlying the jet quenching phenomenon. 

QCD evolution of jets in the quark-gluon plasma 

Svend Domdey 

Based on our recent paper II170L I briefly discuss the effect of scattering for jet quenching of gluonic 
fragmentation functions in the quark-gluon plasma. In contrast to the common discussion Ill631ll6llll71l 
which treats medium-induced gluon emission as the main mechanism, we find that energy loss from 
scattering will also contribute significantly. The basic idea given in the talk is as follows: The DGLAP 
equation is build from the probability for splitting due to change of virtuality. The time scale for this pro- 
cess IS split 

= E /Q^. We extend this equation to contain also scatterings with gluons from the medium 
by constructing a "scattering probability" from the ratio of the time scale for splitting and the scattering 
mean free path igcat = 1/ {'^^g)^ i-C- the number of scatterings which can happen between two splittings. 
Similar to splitting the medium-modified evolution equation contains a gain and a loss term. 

We discuss two slightly different approaches: One is more suitable for analytical calculations, the 
other one for implementation in a Monte Carlo study. Both methods give very similar results. Even with 
a perturbative cross section, we find a suppression of the fragmentation function in medium relative to 
vacuum which is increasing with energy fraction x. For Q =100 GeV and a plasma temperature of T = 
500 MeV (which may be relevant for LHC) we find a suppression of 10% at a: ~ 0.1, but a suppression 
of factor 2 and more for large x. For RHIC conditions, this suppression is more moderate. 

Jet reconstruction performance in the ALICE experiment 

Magali Estienne 

Under extreme conditions of pressure or temperature, nuclear matter is expected to undergo a phase 
transition to a deconfined medium, the quark gluon plasma (QGP). Such an atypical state of matter can 
be defined as a very high density medium in which free partons strongly interact. Experimentally, it is 
possible to re-create the conditions necessary to its formation by colliding heavy ions at ultra-relativistic 
velocity. It is part of the experimental program which will be performed at the Large Hadron Collider at 
CERN. ALICE, one of the four LHC experiments, is dedicated to the study of the QGP properties. 

In QCD, jets are defined as cascades of consecutive emissions of partons initiated by partons from an 
initial hard scattering. This process is followed by the fragmentation of partons into clusters of particles 
emitted in some opposite coUimated directions. In heavy ion collisions, the scene of parton fragmenta- 
tion is modified from vacuum to a dense medium. Before they fragment, partons are expected to lose 
energy through collisional energy loss and medium induced gluon radiation, dominant process in a QGP, 
also called "jet quenching" [171]. Direct consequences of this quenching effect are a modification of 
the jet shape as well as of the parton fragmentation function associated II168I . The ALICE experiment 
is planning to measure these modifications thanks to its excellent tracking capabilities down to momenta 
of order of Agcz) and up to 100 GeV/c and thanks to its large PID capabilities II172II . 
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Different types of jet finding algoritfims (cone and kj) are under study in ALICE for botfi p-p 
and Pb-Pb collisions II1731I1741 . Because of the large underlying background in a jet cone of radius 
R = \J Arj^ + A(|)2 = 1 as well as its large fluctuations in the LHC heavy ion environment, some studies 
have been performed in order, first, to reduce their contributions in each event, and then to subtract the 
remaining contributions event-by-event in the process of jet finding. For the former, HIJING simulations 
have shown that the application of a pr cut of 1 or 2 GeV/c on charged particles and the reduction of the 
jet cone size to R = 0.3-0.4 have to be applied. For the later, three methods of subtraction (statistic, cone 
and ratio) have been tested and compared II1721I175I . With the iterative jet cone finder based on the UAl 
cone algorithm and a full detector simulation, a mean cone energy of (£") = 45 GeV and a resolution 
of RMS/ {E) ~40% are obtained using only charged particles in the reconstruction of 100 GeV mono- 
energetic jets in p-p collisions. The inclusion of neutral particles measured in the ALICE calorime- 
ter HI 761 allows one to increase (E) to 70 GeV while improving the resolution up to ~ 30% (for a jet 
cone radius of R = 0.4). A deeper study of background subtraction in Pb-Pb collisions has shown that the 
jet reconstruction algorithms bias the mean background energy contained in the jet cone toward values 
higher than the ones found in the same area outside the jet. This bias leads to an over estimation of the 
reconstructed energy II177II . We obtain a preliminary (£) « 75 GeV and a resolution of 30-35% for the 
case charged -i- neutral particles for 100 GeV mono-energetic jets in minimum-bias Pb-Pb collisions. 

The full jet spectrum has also been simulated and studied in p-p and Pb-Pb collisions. Due to the 
steeply falling jet production spectrum and the background and signal fluctuations of the reconstructed jet 
energy, we observe a smearing of the reconstructed spectra in p-p. This effect is increased in Pb-Pb due 
to the background fluctuations and its bias in the jet finding process. Jets clearly present advantages for 
the study of jet quenching at the LHC compared to the studies of leading particles which are trigger- and 
surface-biased. With its reconstruction capabilities, ALICE can reconstruct jets with enough resolution 
for Pb-Pb studies even though the observed smearing of the reconstructed spectrum, due to the steeply 
falling shape, will have to be well controlled both in elementary and Pb-Pb collisions. 

Parton fragmentation studies in ATLAS 

Jin Dolejsi, for ATLAS collaboration 

ATLAS installation and commissioning proceeds towards ability to record first pp collisions during 
the year 2008. The ATLAS heavy ion working group has examined the detector performance in condi- 
tions of Pb-Pb collision at 5.5 TeV per nucleon, based on full simulations of the detector response, and 
continues in detailed studies II178II179 1. The focus of this talk is on the jet studies. A dedicated jet finder 
package has been developed and integrated into the standard ATLAS simulation framework. It contains 
the background subtraction, where the background is estimated over the region outside jet candidate, 
and iteratively uses the jet cone algorithm on the calibrated energy depositions in individual cells of the 
ATLAS calorimeters. This package has been tuned on PYTHIA jets embedded into HIJING Pb-Pb colli- 
sion with a hard-cut of high px processes. The simulations show a reconstruction efficiency of the cone 
algorithm steeply rising from about 50% at Ej = 60 GeV up to 90% above 90 GeV. The fake fraction is 
about 0.04 at Ej = 60 GeV and decreases with higher Ej. These results were obtained for events with 
moderately high dN/dr\ = 2700. The axis of the jet can be measured with an angular resolution of about 
0.05 (both in (|) and r\) for jets with Ej = 60 GeV and improves for higher Ej's. The energy resolution is 
about 30% at 60 GeV and is flat over the accessible r\ range. 

The precise position resolution enables a reliable extraction of the jet shapes. The distributions of 7V 
(particle momentum perpendicular to the jet axis) as well as jet fragmentation functions D(z) have also 
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been studied using the tracks reconstructed in the ATLAS Inner Detector. The tracking performance was 
found to be sufficiently good for these distributions. The agreement between the true and reconstructed 
distributions is good, providing that the tracks from the underlying event aie subtracted. The fast-ZTj- 
algorithm [174] is studied as one of the possible alternatives. Its efficiency and energy resolution is 
comparable to that of the cone algorithm and more detailed studies continue. 

Extracting fragmentation functions with y-tagged jet events in Pb-Pb with the CMS experiment 

Christ of Roland, for the CMS Collaboration 

The energy loss of fast partons traversing the strongly interacting matter produced in high-energy 
nuclear collisions is one of the most interesting observables to probe the nature of the produced medium. 
The collisional and radiative energy loss of the partons will modify their fragmentation functions de- 
pending on the path length of the partons in the medium and the medium density. Pb-Pb collisions at 
^/s^ = 5.5 TeV at the LHC will open access to the first detailed measurements of the in-medium modi- 
fications of fragmentation functions of parton initiated jets. The fragmentation functions can be studied 
using the y-jet channel [1 80,1811 . where the initial transverse energy of the fragmenting parton can be de- 
termined from the photon transverse energy, since the photon does not strongly interact with the medium. 

The CMS experiment features large acceptance electromagnetic and hadronic calorimeters, and a 
high precision tracking system well suited to nucleus-nucleus collisions at the LHC [182]. To extract 
fragmentation functions of parton initiated jets in y-jet events, isolated photons reconstructed in the elec- 
tromagnetic calorimeter | < 2.0) are selected and correlated with back-to-back reconstructed calorime- 
ter jets (|ri| < 2.0) in Pb-Pb events. For the selected photon jet pairs, charged particles reconstructed in 
the tracking system (|ri| < 2.5) that lie within a 0.5 radius cone in r\-(^ around the reconstructed jet axis 
are selected. The fragmentation function is constructed correlating the transverse energy of the photons, 
with the transverse momentum of the tracks in the jet cone. 

This study based on a full GEANT simulation of a data set of y-jet events, containing a y with 
Ej > 70 GeV, embedded in central Pb-Pb events corresponding to a full year of data taking with an 
integrated luminosity of 0.5 nb^ ' . The capability to measure jet fragmentation functions in the y-jet 
channel with high statistics is demonstrated. The systematic uncertainties are estimated to be on the 
order of 30% (10%) for quenched events containing a photon with Ej > 70 GeV (100 GeV) mostly due 
to a bias induced by the finite efficiency of the jet reconstruction. Overall the fragmentation function 
can be reconstructed with high significance for 0.2 < ^ < 5 for both, Ej > 70 GeV and > 100 GeV. 
Reconstructed y-jet events can be used to study the dependence of high-pr parton fragmentation on the 
medium and will provide sensitivity to the foreseeable changes in the fragmentation functions relative 
to vacuum fragmentation. This measurement will allow a quantitative test of proposed mechanisms for 
medium-induced parton energy loss, testing fundamental properties of high-density QCD. 

Issues of fragmentation function and medium effects in single inclusive production and two particle 
correlations in p-p and A-A collisions 

Michael Tannenbaum 

Hard-scattering of point-like constituents (partons) in p-p collisions was discovered at the CERN- 
ISR in 1972 in measurements of inclusive single or pairs of hadrons with large transverse momen- 
tum (pr). Due to the steeply falling power-law pj spectrum of the hard-scattered partons, the inclusive 
single particle (e.g. 7i°) pj, spectrum from jet fragmentation of a parton with pj, is dominated by trigger 
fragments with large {zt) ~ 0.7 — 0.8, where Zt = Pt, / Pt, is the fragmentation variable. It was generally 
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assumed, following Feynman, Field and Fox II183I . as shown by data from the ISR experiments, that the 
distribution of away side hadrons from a single particle trigger [with /JrJ, corrected for {zt), would be 
the same as that from a jet-trigger and follow the same fragmentation function measured in e^e^ or DIS. 
PHENIX Ill84[|185ll attempted to measure the fragmentation function from the away side xe ~ pr^ / Pt, 
distribution of charged particles triggered by a in p-p collisions and showed by explicit calculation 
that the xe distribution is actually quite insensitive to the fragmentation function. A new formula for the 
distribution of an associated away-side particle with transverse momentum pj^, which is presumed to 
be a fragment of an away -jet with ^r,, with exponential fragmentation function D(z) = Be^^"-, triggered 
by a particle with transverse momentum pT,, presumably from a trigger-side jet with pT,, invariant cross 
section. EcPoldp^ =Apj". was given HI 8411 : dPp^^/dxE\pj^ ^ (i+Wit)" ' This formula relates jc^, the 
ratio of the transverse momenta of the measured particles, to Xf, = pT^/pr,, the ratio of the transverse mo- 
menta of the away-side to trigger-side jets, where (m) is the mean multiplicity of particles in the away jet. 

Many analyses of the away-jet pj^ distributions in Au-Au collisions are available; but these tend 
to describe the effect of the medium with the variable Iaa{xe), the ratio of the xe distribution in A- A 
collisions to that in p-p collisions, which typically shows an enhancement at low values of xe and a sup- 
pression at higher values of xe- Such behavior could be explained as a decrease in x/, in A-A collisions 
due to energy loss of the away jet in the medium. Fits of the above formula to the available data were 
presented to establish whether: a) the away-jets simply lose energy; b) some of the away-jets lose energy, 
others punch-through without losing energy; etc. 

Measurements of direct single y production in p-p and Au-Au are also discussed. Fragmentation 
photons in p-p collisions, if they exist, can be eliminated by isolation cuts, although, at the ISR and 
at RHIC, measurements of correlations to a direct single y show very little activity, if any, on the same 
side. In contrast to triggers with 71'' which are the fragments of jets, away-side correlations for direct 
single y's which are direct participants in 2-2 hard scattering should measure the fragmentation function 
of the away side jet, modulo kj smearing. This fragmentation function can be plotted conventionally as 
a function of xe or as a function of ^ = ln(l/x£) as suggested for jets at the LHC [il68ll . 

Fragmentation function from direct-photon associated yields at RHIC 

Jan Rak 

The interest of heavy-ion community in the fragmentation function (FF) is manifold. The detailed 
knowledge of the FF modification in heavy-ion collisions is important for constraints of various jet- 
quenching models [1681. Another region of interest is related to the 777 -broadening predicted to accom- 
pany the parton interaction with excited nuclear medium 1 186 1. One way of evaluating the -broadening 
is to measure the parton intrinsic momentum k-^ [184]. It has been shown that in order to extract the value 

of \J{j^ the knowledge of the average trigger particle momentum fraction (zt) and thus FF is vital (see 
Eq. (22) in II184II ). It has been also demonstrated that, despite a common belief, an analysis of charged 
particle distributions associated with the high-/??- trigger hadron cannot provide sufficient constraint for 
the experimental determination of FF. The reason for that luck of sensitivity of high-pj- trigger hadron 
associated spectra to the shape of FF lies in the fact that when events with different associated trans- 
fer momenta are selected the trigger particle momentum fraction (zt) changes event though the trigger 
particle momentum is fixed. In this talk I have presented an analysis of direct-photon associated spec- 
tra which should be free from the above mention problem. There is, however, an effect of ^T-smearing 
which causes the imbalance between the photon-quark momenta. I have presented a method of unfolding 
the /cj-smearing which allows to recover the shape of FF from photon associated spectra. 
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Systematics of complete fragment distributions from nuclear collisions 

Thomas A. Trainor 

The study of fragmentation in RHIC nuclear collisions addresses the central question how formation 
and evolution of a QCD medium in heavy ion collisions may affect parton scattering and fragmentation. 
Whereas previous emphasis has been placed on larger parton energy scales described by perturbative 
methods, we seek a complete description of all fragmentation processes in nuclear collisions. Fragmen- 
tation in elementary collisions provides a reference system. To better determine the energy scale depen- 
dence of fragmentation functions from e+e^, the FFs were replotted on a normalized rapidity variable, 
resulting in a compact form precisely represented by the beta distribution, its two parameters varying 
slowly and simply with parton energy scale Q [187|. The new parameterization extrapolates fragmenta- 
tion functions to small Q and describes the fragmentation process down to small transverse momentum 
Pt in nuclear collisions at RHIC. Analysis of transverse momentum pj spectra to pj = 6 GeV/c for 
10 multiplicity classes of p-p collisions at ^/s = 200 GeV revealed that the spectrum shape can be de- 
composed into a part with amplitude proportional to multiplicity and described by a Levy distribution 
on transverse mass mj, and a part with amplitude proportional to multiplicity squared and described 
by a Gaussian distribution on transverse rapidity jt 111881 . The functional forms of the two parts are 
nearly independent of event multiplicity, and they can be identified with the soft and hard components 
of a two-component model of p-p collisions: longitudinal participant nucleon fragmentation and trans- 
verse scattered parton fragmentation. Interpretation of the hard component as parton fragmentation is 
supported by correlation studies l|189[|190i . 

A similar two-component analysis of spectra to = 12 GeV/c for identified pions and protons from 
200 GeV Au-Au collisions applied to Au-Au centrality dependence revealed that the soft-component 
reference is again a Levy distribution on transverse mass mj, and the hard-component reference is again 
a Gaussian on yj, but with added exponential {pj power-law) tail [191]. Deviations of data from the 
two-component reference are described by hard-component ratio r^A which generalizes nuclear modi- 
fication factor Raa- Centrality evolution of pion and proton spectra is dominated by changes in parton 
fragmentation. The structure of vaa suggests that parton energy loss produces a negative boost Ajj- of a 
large fraction (but not all) of the minimum-bias fragment distribution (hard component), and that lower- 
energy partons suffer relatively less energy loss, possibly due to color screening. 

A recent analysis of azimuth correlations based on methods in II192II193TI has isolated azimuth 
quadrupole components from published V2{pt) data (called elliptic flow) as spectra on transverse ra- 
pidity yj for identified pions, kaons and lambdas from minimum-bias Au-Au collisions at 200 GeV. The 
form of the spectra reveals that the three hadron species are emitted from a common source with boost 
^yto ~ 0.6. The quadrupole spectra have a Levy form on similar to the soft components of single- 
particle spectra, but with significantly reduced (~ 0.7 x) slope parameter T . Comparison of quadrupole 
spectra with single-particle spectra suggests that the quadrupole component comprises a small fraction 
(< 5%) of the total hadron yield, contradicting the hydrodynamic picture of a thermalized, flowing bulk 
medium. The form of V2{pt) is, within a constant factor, the product of p'^ {pj in the boost frame) times 
the ratio of the quadrupole spectrum to the single-particle spectrum. The form of the spectrum ratio then 
implies that above 0.5 GeV/c, V2{pt) is dominated by the hard component of the single-particle spectrum 
(minijets). It is thus unlikely that so-called constituent-quark scaling attributed to V2 is relevant to soft 
hadron production mechanisms (e.g., chemical freezeout). These results suggest instead that "elliptic 
flow" may result from fragmentation of gluonic quadrupole radiation, not hydrodynamic evolution. 
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